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ESTABLISHMENT  OF  A  ROTOR  MODEL  BASIS 


R.  E.  McFarland 

Ames  Research  Center 
Aeromechanics  Laboratory 
and 

AVRADCOM  Research  and  Technology  Laboratories 


I.  SUMMARY 


Rotating  blade-element  mathematical  models  for  use  in  discrete,  man-in-the-loop 
rotorcraft  simulation  have  considerable  engineering  value,  but  they  are  computation¬ 
ally  expensive.  A  comprehensive  blade-element  model  of  the  RSRA  rotor  system  is 
examined  here  to  provide  accurate  baseline  data  for  the  analysis  of  modifications  to 
the  present  model.  It  is  important  that  models  be  developed  for  real-time  simulation 
which  retain  the  important  features  of  the  original,  while  significantly  reducing  the 
computational  expense. 

Lift  and  drag  data  are  presented  for  each  of  five  blades  at  aircraft  speeds 
up  to  300  knots.  Total  main  rotor  forces  and  moments  are  also  provided  over  this 
speed  range.  This  body  of  data  is  the  basis  for  comparison  with  data  from  alternate 
rotor  models. 

Obtaining  these  data  required  developments  in  both  rotorcraft  modeling  and 
computer  sciences.  In  the  first  discipline  it  was  determined  that  the  model  for  the 
blade  element  angle  of  attack  should  be  modified.  This  modification  was  used  in  the 
generation  of  the  basis  data.  Also,  the  tip-loss  model  was  found  to  require  modifi¬ 
cation  when  a  large  number  of  segmental  elements  is  used,  as  in  the  case  of  a  basis 
generation.  An  exploration  of  this  phenomenon  leads  to  an  inverted  tip-loss 
definition. 

The  computer-science  impact  is  evaluated  in  the  discussion  of  elemental  spacing 
criteria;  this  evaluation  influences  the  method  of  integrating  the  elemental  forces 
and  moments  along  the  radius  of  a  rotor  blade.  It  is  shown  that  locations  within 
segments  at  which  the  lift  and  drag  evaluations  are  made  are  critical  to  accuracy 
when  the  quadrature  technique  is  used,  and  that  inconsistent  forces  and  moments  may 
result.  This  aspect  is  examined,  and  independent  of  the  spacing  algorithm,  consistent 
summation  coefficients  are  developed  using  the  assumption  of  linear  force  profiles 
between  observation  points.  Another  computer-science-related  idea  which  is  developed 
involves  the  technique  of  performing  the  force  and  moment  evaluations  at  fixed 
azimuthal  positions  for  ease  in  data  comparison.  Coupled  with  this  idea  is  the 
observation  that  the  aerodynamic  and  kinematic  portions  of  the  dynamics  of  a  blade 
may  then  be  independently  analyzed.  The  idea  of  "constrained"  dynamics  allows 
aerodynamic  parametric  changes  to  be  evaluated  by  essentially  holding  the  kinematic 
portions  to  a  prescribed  condition.  The  result  is  a  considerable  reduction  in 
computation  required  for  data  collection,  and  a  clear  establishment  of  cause  and 
effect  for  analysis. 

A  very  accurate  basis  for  model  comparison,  which  is  appropriate  for  explora¬ 
tions  into  real-time  model  efficiencies,  has  emerged  from  this  study.  Many  of  the 
corrections  and  improvements  developed  are  recommended  for  immediate  applicability 
to  a  simulation  model  already  being  used.  Although  the  spacing  characteristics  of 


the  equal  annuli  algorithm  are  preserved,  alternate  techniques  are  recommended  for 
the  use  of  summation  coefficients,  and  the  computation  of  the  tip-loss  phenomenon. 


II.  INTRODUCTION 


The  RSRA  simulation  mathematical  model  was  developed  by  Sikorsky  Aircraft  for 
performance  and  handling  qualities  evaluations;  it  is  also  intended  "for  pilot  train¬ 
ing,  preflight  of  test  programs,  and  the  evaluation  of  promising  concepts  before 
their  implementation  on  the  flight  vehicle"  (Houck  et  al.,  1977).  The  author  is 
indebted  to  Houck  et  al.,  for  the  extensive  use  of  their  data  in  this  study.  It  is 
assumed  that  the  reader  has  some  familiarity  with  blade-element  mathematical  model¬ 
ing,  and  in  particular  with  the  main  rotor  model  described  by  Houck  et  al.,  which  is 
herein  referred  to  as  the  "extant  model."  As  supplied  to  the  Langley  and  Ames 
Research  Centers,  the  model  represents  a  computational  requirement  that  can  only  be 
synchronized  with  real  time  by  the  use  of  rather  sophisticated  computers.  Because 
of  the  applicability  of  the  blade-element  model,  NASA  has  an  interest  in  its  wider 
utility  (Mackie  and  Alderete,  1977),  for  instance,  on  such  modest  computers  as  the 
Xerox  Sigma  8.  The  equations  and  data  contained  in  the  real-time  model  represent  a 
proprietary  derivation  from  Sikorsky  Aircraft's  General  Helicopter  (GENHEL)  simula¬ 
tion  model.  Only  the  radial  dimension  is  investigated  here. 

The  objective  of  this  study  is  to  establish  a  rotor-system  data  base  of  high 
fidelity,  using  relatively  arbitrary  inputs  that  may  be  used  in  establishing  the 
accuracy  of  alternate  mathematical  models.  It  is  assumed  that  these  alternate  math¬ 
ematical  models  will  use  the  same  physical  principles  as  the  extant  rotating,  blade- 
element  model,  but  they  will  contain  various  modifications  in  order  to  accrue  real¬ 
time  computational  benefits.  Because  of  the  magnitude  of  the  computational  problem 
for  blade-element  rotor  models,  their  accuracy  is  highly  correlated  with  their  execu¬ 
tion  speed.  It  Is  maintained  that  an  approximate  model  may  be  developed  that  will 
compare  favorably  with  the  complete  model  by  capitalizing  upon  speed  advantages.  If 
these  advantages  approach  the  factor  of  two,  the  model  will  then  be  capable  of  real¬ 
time  synchronization  using  a  modest  digital  computer  such  as  the  Xerox  Sigma  8. 

Hence,  real-time  constraints  are  the  motivation  for  this  study,  although  only  accu¬ 
racy  is  addressed  here. 

In  order  to  establish  a  basis,  or  data  base  of  high  fidelity  for  comparison  pur¬ 
poses,  it  must  be  assumed  that  the  selected  discretized  model  collapses  to  the  proper 
continuum  when  the  number  of  discrete  elements  is  large  (but  considering  computer 
word  length,  not  too  large).  Unfortunately,  the  extant  RSRA  simulation  model  (Houck 
et  al. ,  1977)  converges  very  slowly  with  segment  count,  and  this  behavior  prompted 
an  investigation  of  the  radial  spacing  convergence  properties.  Both  a  low-  and  high¬ 
speed  flight  condition  were  selected,  and  are  presented  in  figures  1  and  2.  Curves 
tagged  "extant"  reveal  the  convergence  difficulties.  The  rotor-system  thrust  differ¬ 
ences  as  functions  of  the  total  number  of  evaluation  points  are  slowly  convergent  for 
a  large  number  of  points,  and  radically  divergent  for  a  small  number  of  points. 

Hence,  in  this  document  we  shall  go  somewhat  beyond  the  mere  creation  of  a  data  base 
and  develop  certain  corrections  and  improvements  that  are  applicable  to  the  extant 
real-time  simulation  model. 

A  basis  is  created  which  meets  the  test  of  invariancy  with  perturbations  in  the 
number  of  radial  evaluations.  This  corrected  and  improved  model  is  not  for  real-time 
use,  but  for  establishing  a  standard,  for  comparison  purposes,  with  20  radial  evalua¬ 
tion  points.  This  "basis"  produces  the  rotor-system  thrust  differences  tagged  "basis" 
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in  figures  1  and  2.  The  use  of  20  evaluation  points  clearly  eliminates  discretiza¬ 
tion  phenomena  and  algorithmic  distortion  from  this  model. 


In  contrast  to  the  extant  real-time  model  (Houck  et  al.,  1977),  an  alternate 
model  is  developed  which  has  better  convergence  features.  This  model  is  also  appro¬ 
priate  for  real-time  simulation,  and  delivers  the  thrust  differences  tagged  "alter¬ 
nate"  in  figures  1  and  2.  Since  real-time  constraints  usually  dictate  a  small 
number  of  evaluation  points,  this  alternate  model  is  recommended  for  simulation  work. 

In  the  process  of  creating  a  data-acquisition  algorithm  for  this  analysis,  a 
technique  was  discovered  for  literally  taking  a  snapshot  of  dynamic  va~iables;  this 
technique  is  useful  in  isolating  the  small  differences  normally  associated  with  vari¬ 
ous  discrete  algorithms.  This  feature  has  considerably  reduced  the  computational 
expense  of  this  and  other  studies;  and  it  is  also  useful  for  the  problem  of  rotor- 
craft  trimming. 

A  technique  is  developed  for  handling  the  tip-loss  phenomenon  by  introducing 
the  concept  of  a  pseudotip  for  lifting  purposes.  This  technique  reduces  computa¬ 
tional  errors  and  simplifies  the  generalization  of  the  integration  process  to  a  sum¬ 
mation  process. 

With  respect  to  a  force  profile  along  a  blade,  the  combination  of  summation  and 
spacing  algorithms  are  investigated  and  found  to  be  less  than  optimum.  Improvement 
to  the  summation  algorithm  motivates  the  development  of  a  slightly  different  spacing 
algorithm,  used  primarily  for  illustrative  purposes;  it  is  identified  as  equal- 
annuli  midpoint  spacing  and  is  shown  to  better  accommodate  the  quadrature  technique. 
Without  influencing  computational  expense,  summation  coefficients  are  developed  for 
improvements  in  accuracy,  and  for  consistency  between  forces  and  moments,  regardless 
of  the  spacing  algorithm.  These  coefficients  feature  a  linear  assumption  for  elemen¬ 
tal  forces  between  evaluation  radii  and  extrapolation  at  endpoints. 

A  five-point  model  that  uses  the  corrections  and  improvements  developed  here 
becomes  so  accurate  that  it  is  difficult  to  justify  the  use  of  more  segments  in 
real-time  simulation  work. 


III.  DATA  ACQUISITION 


An  extensive  mathematical  model  of  the  RSRA  (Houck  et  al.,  1977)  has  been  pro¬ 
grammed  at  NASA/ Ames  Research  Center  in  both  a  single  and  a  dual  computer  configura¬ 
tion  (Mackie  and  Alderete,  1977).  Although  satisfactory  real-time  synchronization 
with  this  software  has  not  been  achieved  using  modest  digital  computers,  the  single 
computer  formulation  is  nonetheless  readily  available  for  analyses  using  data  gen¬ 
erated  in  a  batch  (nonreal-time)  mode.  When  an  appropriate  driving  routine  is  used, 
the  subroutine  identified  as  ROTOR  may  be  exercised  as  a  separate  entity  for  investi¬ 
gations  into  rotor  modeling  techniques. 

In  order  to  create  a  basis  with  which  to  compare  alternate  model  structures, 
extensive  data  were  obtained  from  ROTOR.  Using  20  segments  per  blade,  and  with  modi¬ 
fications  that  will  be  developed  herein,  ROTOR  was  run  with  an  azimuth  advance  angle 
of  9  ;  the  (pseudo)  cycle  time  was  varied  from  7.094  to  9.851  msec  in  order  to  accom¬ 
modate  the  scheduled  rpm  constraints  (Mach  1  avoidance)  as  shown  in  table  1. 


Data  acquisition  occurred  at  multiples  of  eight  times  the  9°  azimuth  advance 
increment,  that  is,  every  72°,  in  order  to  avoid  data  interpolation. 

For  an  arbitrary  input  control  configuration  within  the  capacity  of  a  rotor 
system,  after  blade  transients  vanish  after  a  few  revolutions,  the  tracking  blade 
phenomenon  is  observed  in  this  model  such  that  at  any  selected  azimuth  station  the 
dynamics  of  a  given  blade  are  replicated  by  the  other  blades  as  they  arrive  at  that 
station.  In  this  condition  of  equilibrium,  blade-element  rotational  kinematics  at 
these  orientations  may  be  described  with  Fourier  coefficients.  For  the  RSRA  model 
at  all  flight  regimes  this  phenomenon  has  been  observed  to  occur  in  something  less 
than  ten  complete  rotor  revolutions,  or  400  cycles  through  the  9°-increment  logic. 
The  orientation  given  in  figure  3  has  been  selected  as  a  standard  for  observation 
purposes,  where  the  "master  blade"  is  pointing  aft. 


The  problem  of  obtaining  a  control  profile  that  is  within  the  capacity  of  the 
rotor  system  was  treated  with  some  seriousness,  although  it  is  not  necessarily  ger- 
main  to  this  study.  All  that  is  actually  required  in  this  analysis  is  that  all  the 
degrees  of  freedom  be  exercised,  and  this  requirement  turns  out  to  be  a  foregone 
conclusion  with  an  arbitrary  control  profile.  Nonetheless,  a  voluminous,  unpublished 
document1  was  consulted  for  data  over  the  entire  flight  envelope,  and  a  control  pro¬ 
file  was  selected  that  produced  a  reasonably  smooth  transition  from  point  to  point  in 
velocity  space  (10-knot  increments).  This  profile,  shown  in  figure  4,  is  sufficient 
to  acquire  a  large  data  base  from  the  ROTOR  routine,  as  modified.  Since  the  rotor  as 
a  system  has  been  isolated  in  this  study,  it  is  not  necessary  to  consider  the  viabil¬ 
ity  of  the  selected  control  profile  with  respect  to  body-loop  closure.  These  controls 
certainly  do  not  constitute  trim  positions,  and  the  probability  or  even  possibility 
that  they  will  actually  occur  "in  flight"  is  questionable. 


It  is  shown  herein  how  tip-loss  effects  may  be  independently  considered  by  util¬ 
ization  of  the  concept  of  a  pseudotip  for  lifting  purposes,  and  how  this  concept 
eliminates  discontinuous  end  effects  from  the  force  profiles  as  a  function  of  radius. 
The  integration  (summation)  of  these  smoother  "force  profiles," 


fp  =  dFp/dr  ' 
f*p  =  dFj/dr 
fR  -  dFR/dr  , 


(1) 


is  shown  still  to  present  an  array  of  discretization  problems,  which  are  compounded 
when  the  aerodynamic  moments  are  also  considered. 


With  two  fixed  locations  along  the  blade,  one  being  the  cuff  (beginning  of  the 
effective  blade)  and  the  other  being  the  pseudotip  (which  is  shown  to  be  the  effec¬ 
tive  tip  for  lifting  purposes),  the  blade-force  profiles  are  presented  by  using 
20  evaluations  per  blade,  and  these  data  are  given  every  10  knots  from  0  to  300  knots; 
lift  profiles  are  presented  in  figure  5,  and  drag  profiles  are  presented  in  figure  6. 
The  quantity  fR  is  a  minor  term  and  is  not  displayed  here. 

All  the  data  in  this  document  have  been  created  after  the  implementation  of  the 
angle-of-attack  correction  given  in  the  next  section. 


Robert  A.  Monteleone,  Systems  Requirements  Handbook  for  the  Rotor  Systems 
Research  Aircraft,  Sikorsky  Report  SER-72039,  March  1977.  Also  available  from  NASA 
management,  Contract  NAS1-13000. 
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Angle-of-Attack  Error 


In  the  reference  formulation  (Houck  et  al. ,  1977)  ay  is  a  four-quadrant  angle 
(180°  <  a  £  180°)  that  may  be  given  in  terms  of  the  dimensional  velocity  variables 
(vp,  vT,  vR)  by 


ay  =  tan" 


(vy  tan  6^  +  vp)|cos  y| 
VT  ~  VP  tan  ®A  cos2  Y 


When  the  denominator  of  this  expression  is  positive,  | ocy |  £  90°  and  the  positive 
table  lookup  argument  (<*TRANS  for  lift)  is  obtained  by  determining  whether  or  not 
the  inplane  angle  of  attack  given  by  | ay  cos  y|  is  within  a  small  band  about  the 
positive  abscissa: 


“TRANS  = 


When  the  denominator  of  equation  (2)  is  negative,  | oty |  >  90°  and  reverse  flow 
occurs.  The  value  of  <*TRANS  for  lift  is  then  determined  by  examining  a  similar 
small  band  about  the  negative  abscissa: 


aY  cos  y | 

( | ay  cos 

Yl 

^  13.5°, 

1  “Y  1 

1  A 

vO 

o 

o 

v-/ 

“yI 

(  j  Oty  COS 

y| 

>  13.5°, 

1  “Y 1 

1  A 

VO 

o 

o 

“TRANS  = 


180°  -  (180°  -  |  ay  |  )  |  cos  y|  [(180°  -  |aY|)|cos  y|  <  8°,  [oty  1  >  90°] 
| “y |  [(180°-  |aY|)|cos  y|  >  8°,  | oty |  >  90°] 


Section  A-20(c)  of  the  paper  by  Houck  et  al.,  1977  (the  first  line  of  eq.  (9)  below) 
is  therefore  assumed  to  be  in  error,  and  should  appear  as 

&Y 

“TRANS  “  “y I  cos  y |  +  -  180° (1  -  | cos  y|) 

I  “Y I 

=  180°  -  (180°  -  |ay | ) [cos  yj  (5) 

which  is  the  same  as  the  first  line  of  equation  (4).  Since  the  angle  of  attack  is 
also  used  to  compute  the  switch  point  between  univariant  and  multivariant  lift  and 
drag  functions,  this  section  of  the  function-generation  logic  is  clarified  in 
figure  7. 

This  angle-of-attack  error  in  reverse  flow  probably  has  little  effect  on  previ¬ 
ous  simulation  results,  although  it  has  caused  considerable  delay  in  this  analysis. 
Its  influence  is  examined  below. 

If  the  equal-annuli  spacing  algorithm  of  the  extant  model  is  used,  the  initial 
evaluation  (inboard)  radius  is  given  by  (see  section  VII) 

/r2  -  r2 

1/  MR  A  .  o 

ri  *  V - 2N -  +  rA  <6) 
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which  for  five  segments  is  equal  to  11.56  ft.  The  interior  portion  of  the  effective 
aerodynamic  blade,  which  is  therefore  treated  as  lacking  in  any  important  functional¬ 
ity,  is  given  by  the  proportion 


P  = 


rl  -  rA 
BMRrMR  "  rA 


(7) 


or  21.6%  for  a  five-segment  simulation.  For  the  simplified  relationship  of  tangential 
velocity. 


vT  =  «r  +  1.69(vknots)sin  ip 


(8) 


at  the  worst-case  azimuth  of  270°,  and  setting  the  tangential  velocity  equal  to  zero 
permits  the  calculation  of  the  radius  at  which  stall  occurs;  this  calculation  is  pre¬ 
sented  in  figure  8.  If  the  computed  values  of  rpm  of  table  1  are  used,  it  is  seen 
that  the  stall  phenomenon  begins  to  occur  on  the  effective  aerodynamic  surface  at 
about  85  knots  of  vehicle  velocity.  It  is  also  seen  that  until  about  155  knots  this 
phenomenon  is  ignored  in  the  extant  model  since  it  occurs  interior  to  the  initial 
evaluation  radius  (five  segments).  Hence,  even  without  the  angle-of-attack  correc¬ 
tion  given  here,  the  extant  simulation  model  performs  as  it  was  intended  in  velocity 
regions  less  than  155  knots.  Beyond  this  velocity,  or  if  more  than  five  segments  are 
used  at  lower  velocities ,  the  errors  are  dramatic. 

As  an  example  of  the  effect  of  this  error,  let  us  consider  the  extant  (erroneous) 
formulation,  which  is  independent  of  the  sign  of  ay  whenever  the  absolute  value  of 
ay  is  greater  than  90° : 


“TRANS 


“Y 


cos  y  + 


180° (1 


I  cos  y|) 


cos  y|«yI  +  180° (1  -  | cos  y| ) 


=  I  — | ay  cos  y|  +  180° (1  -  | cos  y J ) 

The  latter  step  is  performed  when  |ay|  >  90°  because  cos  y  is  negative  when 

v>p  <  0 . 


(9) 


If  ay  =  179°,  which  is  well  within  the  86  band  about  the  negative  abscissa, 
then  “TRANS  should  be  even  closer  to  180*.  Also,  for  cos  y  "  0.9,  where  its  sign 
is  immaterial,  one  would  anticipate  a  value  for  “TRANS  somewhere  within  1°  of  the 
negative  abscissa  (180°),  but  what  actually  occurs  from  equation  (9)  is: 


“TRANS 


- 1  (179°) (0. 9) |  +  180° (1  -  1 0 . 9 1 ) 


143.1° 


(10) 


The  correct  value  from  equation  (5)  is  179.1°  (closer  to  coincidence  with  the 
airstream),  which  produces  a  relatively  small  lift  coefficient  (about  -0.09).  The 
value  143.1°  produces  a  huge  lift  coefficient  (about  -1.1)  and  has  been  observed  to 
cause  spectacular  N-per-revolution  force  and  moment  variations. 


Data  Density 


The  blade  orientation  of  figure  3  and  the  control  profile  of  figure  4  were  used 
in  order  to  obtain  a  reasonable  output  history;  figures  5  and  6  show  the  diverse 
albeit  piecewise  continuous-curve  shapes  of  the  individual-blade  lift-  and  drag-force 
profiles.  Twenty  observations  of  each  lift-  and  drag-profile  curve,  as  shown  in 
these  figures,  appear  to  be  more  than  heuristically  sufficient  for  continuum  emula¬ 
tion,  and  this  fact  is  reinforced  by  invariancy  under  summaLion,  as  is  shown  by  the 
curves  tagged  "basis"  in  figures  1  and  2.  (The  actual  summation  algorithm,  which 
produces  invariant  answers  for  a  reasonable  number  of  segments  N,  will  be  seen  to  be 
different  if  a  minimal  value  of  N  is  required.)  However,  with  the  luxury  of 
20  observations,  the  summation  algorithm  becomes  relatively  Inconsequential. 

The  data  displayed  in  figures  5  and  6  contain  the  angle-of-attack  correction 
discussed  previously.  The  observation  radii  that  have  been  selected  for  this  data 
display  are  uniformly  spaced,  and  although  the  inboard  observation  is  placed  at 
the  beginning  of  the  aerodynamic  surface  of  the  blade,  the  outboard  observation 
radius  will  be  seen  to  have  a  unique  definition  in  a  subsequent  section  of  this 
report. 

It  has  previously  been  ascertained  that  "reducing  blade  segments  does  not  appear 
to  influence  the  solution  to  any  great  extent  (Houck  and  Bowles,  1976).  Thus,  one 
might  question  the  necessity  for  such  a  fine  mesh  of  points  in  this  analysis.  In 
order  to  answer  this  question,  let  us  consider  the  following: 

In  this  study,  under  all  conditions,  only  a  very  small  azimuth  advance  angle  of 
9°  is  used,  thereby  effectively  eliminating  this  degree  of  freedom  as  an  error 
source.  It  has  indeed  been  our  experience  that  "the  worst  single  effect  is  that  of 
increasing  integration  [azimuth]  interval"  (Houck  and  Bowles,  1976).  It  is  also  a 
very  complex  problem,  which  is  under  investigation  but  not  addressed  here.  Only 
improvements  to  the  radial  dimension  are  addressed  here,  and  in  this  dimension,  which 
is  displayed  in  figures  5  and  6,  the  "frequency  content"  is  quite  low,  provided 
that  the  tip-loss  technique  of  section  IV  is  used  (discontinuities  are  avoided). 
Herce,  the  standard  technique  of  reducing  element  size  until  invariant  answers  ensue 
is  quite  proper,  and  figures  1  and  2  indicate  that  the  value  of  20  evaluation  points 
is  also  quite  proper. 


Balanced  and  Constrained  Dynamics 

The  two  conditions,  "balanced  dynamics"  and  "constrained  dynamics,"  are  con¬ 
venient  for  analysis.  In  the  balance  condition  the  flap/lag  moments  are  allowed  to 
operate  as  forcing  functions  for  the  flap/lag  differential  equations  (normal  opera¬ 
tions),  and  hence  produce  new  flapping  and  lagging  states.  In  this  condition,  how¬ 
ever,  the  slight  modification  of  any  parameter  or  algorithm  generally  produces  myriad 
state  changes  and  masks  any  causal  relationship  under  investigation.  For  this  rea¬ 
son,  once  the  tracking-blade  phenomenon  is  observed,  the  applicable  Fourier  coeffi¬ 
cients  are  captured  for  use  in  the  constrained  condition.  Thus,  the  condition  of 
"constrained  dynamics"  means  that  blade  trajectories  are  constrained  to  follow  those 
previously  determined  (by  the  20-segment  model)  so  that  force  and  moment  differences 
may  be  isolated  to  differences  in  algorithms  and  segmentation.  This  process  yields 
for  both  the  balanced  and  constrained  cases  identical  dynamics  at  the  specific 
azimuth  stations  selected  for  observation.  (The  standard  orientation  of  figure  3  is 
always  used  here.)  In  this  orientation,  where  the  master  blade  is  aft  and  others 


follow  sequentially  every  2ir/B  rad  in  the  direction  of  positive  rpm  (B  is  the 
number  of  blades),  flap/lag  states  are  the  same  in  both  the  balanced  and  constrained 
conditions. 

It  is  interesting  to  note  that  in  the  constrained  condition  only  one  computa¬ 
tional  pass  through  the  rotor  equations  is  necessary  in  order  to  determine  control 
variation  responses.  Responses  to  collective  (0CUFF)  an<*  cyclic  controls  (Ajg,  Bis) 
may  be  approximated  by  use  of  partial  derivatives  with  respect  to  low-order  Fourier 
coefficients.  As  a  subject  for  future  research  the  static  feature  of  the 
"constrained-dynamics"  condition  suggests  an  efficient,  general  trimming  algorithm 
for  rotorcraft  since  control-variation  responses  may  then  be  isolated  from  the  sensi¬ 
tive  balance  process.  When  such  a  two-pass  trimming  algorithm  is  performed  the  rela¬ 
tively  long  transients  induced  by  control  variations  (due  to  the  low  damping  in  blade 
responses)  are  constrained  until  approximately  correct  control  signals  for  trim  are 
determined;  the  blades  are  then  released  until  a  balance  is  obtained.  This  process 
is  then  repeated  until  the  constrained  results  become  quiescent.  It  is  believed  that 
this  sequence  is  more  efficient  than  the  alternate  process  of  determining  total 
responses  to  each  control  variation  since  long-period  rotor  transients  interfere  with 
cause  and  effect. 


IV .  TIP  LOSS 


The  tip-loss  technique  is  discussed  and  arguments  are  presented  for  its  modifi¬ 
cation  based  upon  consistency,  applicability,  accuracy,  and  ability  to  remain  rela¬ 
tively  invariant  with  changes  in  the  number  of  segments  N.  Tip  loss  is  presented 
early  in  this  discussion  because  it  results  in  a  complete  redefinition  of  integral 
limits,  which  influences  other  parameters.  The  resultant  concept  permits  the  imple¬ 
mentation  of  tip  loss  as  a  superimposed  quantity  that  may  then  be  ignored  throughout 
the  remainder  of  this  paper. 

Where  r^R  =  31  ft  is  the  total  blade  length,  the  radius  to  the  cuff  from  the 
hub  (rA  =  6.45  ft)  is  given  as  the  sum  of  the  hinge  offset  (e  *  1.05  ft)  and  spar 
(e'  =  5.4  ft): 


rA  =  e  +  e'  =  rMR(C  +  V)  (H) 

In  terms  of  dimensional  segments  Arn  or  nondimens ional  segments  Ayn,  the  sum  over 
all  N  segments  totals  the  effective  blade  distance,  that  is, 

N  N 

£  Arn  =  rMR  £  A?n  "  rMR<X  "  «  "  *’>  *  rMR  “  rA  <12> 

n=  1  n=  1 


or  24.55  ft.  Houck  et  al.  state  that  for  the  outboard  segment  the  lift  coefficient 
is  modified  for  tip  loss  by  use  of  the  relationship 


C 


LY 


f (“TRANS’ 


Mach) 


(13) 


which  is  equivalent  to  an  Nth  segment  radial-distance  scale,  or  interpolation  in 
the  Nth  segment  via  the  proportion 


(14) 


P 


1  - 


*  -  BMR 

AyN 


d  ■  ®mr) 


rMR 

ArN 


where  =  0.97.  Therefore,  the  contribution  of  the  Nth  (outboard)  segment  to  the 

force  summation  process  for  all  terms  that  contain  the  quantity  Cly  (lift  terms)  is 
scaled,  and  from  the  associative  law  this  scaling  is  shown  below  to  be  identical  to 
a  radial-distance  scale. 

Let  us  consider  the  discretized  lift-force-profile  definition  (Houck  et  al., 
1977)  with  its  radial  increment  multiplier  Arn, 

AFpn  -  \  PcvYn((CLYnvTn)/|cos  yni  +  cDYnvpJArn  (15) 


ISgn(aYn)fn  (1  -  n  <  H) 

PSgn(aYn)fn  (n  =  N) 


(16) 


N  N-l 

FP  *  £  ^  Arn  +  £  «n  Arn  +  «N<P  ArN> 
n»  1  n=>  1 

which  is  the  discrete  form  of  the  continuous  integration  process 


Fp 


/•rMR 


h(r)dr  + 


/ 


BMRrMR 


g(r)dr 


(17) 


where  fn  is  a  tabled  function,  and  Sgn(  )  is  plus  or  minus  unity.  For  the  purpose 
of  simplification,  we  may  define 

gn  =  \  PcvYnSgn(aYn)fnvTn/|cos  yn| 
hn  ■  f  PcvYncDYnvPn 

so  that  equations  (15)  and  (16)  may  be  written 

AFPn  “  8n  Arn  +  Arn  d  -  n  <  N> 
and,  including  the  scale  factor  P  in  the  applicable  term, 

AFPN  =  ®N^P  ArN^  ^N  ArN  = 

The  blade  lift  force  is  obtained  by  the  summation  over  all  segments. 


(18) 


(19) 


(20) 


(21) 


Hence,  CLY  scaling  is  equivalent  to  scaling  the  upper  radial  limit  for  terms  involv¬ 
ing  ClY  in  the  force-profile-integration  process.  Even  if  a  more  complex  tip-loss 
formulation  is  assumed,  such  as  linearly  decreasing  ClY  in  some  radial  interval 
within  the  Nth  segment,  the  discrete  force  summation  (with  Bj^  adjustment)  is 
identical.  One  cannot,  therefore,  ascertain  the  motivation  for  CLY  scaling  from 
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the  provided  information  (eq.  (13)),  but  the  effect  is  clearly  the  abrupt  vanishing 
of  cLY"related  terms  at  Bj^r^  as  in  equation  (21).  The  availability  of  the  con¬ 
tinuum  formulation  (eq.  (21))  is  important  when  a  large  number  of  segments  are  con¬ 
sidered,  as  will  be  seen.  The  upper  radial  limit  Bj^r^  will  be  exploited  in  this 
section,  but  it  is  important  to  keep  in  mind  that  both  Clyn  scaling  and  Arjj  scal¬ 
ing  are  identical  techniques  for  the  force  summation  problem,  with  the  exceptions 
that  are  developed  below. 

When  the  number  of  segments  N  is  allowed  to  become  large,  inconsistencies  in 
the  mathematical  algorithms  become  apparent.  One  such  inconsistency  in  the  extant 
model  involves  the  summation  of  moments,  where  the  moment  arm  to  the  proportioned 
segment  is  not  adjusted.  It  is  not  aerodynamically  sound  to  consider  Cly  as  abso¬ 
lutely  constant  in  the  outboard  segment  all  the  way  to  the  tip  r^,  even  if  it  has 
been  conceived  as  proper'1/  scaled  (although  its  actual  characteristic  is  an  open 
question),  so  such  an  adjustment  to  the  moment  arm  is  always  necessary.  Thus,  no 
physical  assumption  exists  where  the  moment-summation  element  in  the  Nth  segment 
reduces  to  the  product  fN  ArNrN  of  the  extant  model. 

In  this  paper  the  discontinuous  cutoff  of  lifting  potential  at  B^r^  is  gen¬ 
eralized  and  assumed  to  be  the  motivation  for  the  formulation  (13).  This  assumption 
permits  the  development  of  a  moment -summation  algorithm  that  is  consistent  with  the 
force-summation  algorithm  without  added  computational  expense.  Also,  this  assumption 
will  lead  to  a  more  general  treatment  of  the  force  profiles  themselves,  which  will 
become  quite  important  when  the  asymmetrical  distribution  properties  of  the  equal- 
annuli  spacing  algorithm  are  considered.  Unless  the  lift  coefficient  is  constant  in 
the  outboard  segment,  the  proportion  operation  (14)  takes  liberties  with  the  mean- 
value  theorem  for  integrals,  and  it  actually  contributes  to  force  as  well  as  moment 
errors. 

Let  us  consider  the  spacing  algorithm  for  the  extant  model  developed  in  sec¬ 
tion  VII.  Using  the  outboard-segment  increment  ArN, 


the  proportion  (14)  goes  negative  when  N  >  16.  This  operation  is  equivalent  to  sub¬ 
tracting  more  than  the  entire  force  contribution  of  the  outboard  segment,  and  in 
moment  space  such  a  subtraction  cannot  be  justified  even  if  the  force  increments  are 
constants.  Indeed,  figure  5  demonstrates  that  there  is  considerable  variation  in 
the  lift-force  profiles  in  this  region,  and  comparison  with  figure  6  shows  that 
typical  lift  values  a r'e  an  order  of  magnitude  larger  than  drag  values.  Since  this 
unusual  subtraction  feature  occurs  just  when  the  quantity  N  is  getting  large 
enough  in  an  accuracy  analysis  to  replicate  continuum  results,  a  large  number  of 
segments  N  and  this  technique  for  tip-loss  compensation  are  inconsistent.  Hence, 
at  least  for  the  establishment  of  a  basis,  this  technique  begs  modification.  It 
should  be  noticed  that  the  outboard-segment  size  Ar^  in  equation  (22)  is  approx¬ 
imately  hyperbolic  with  N  (Taylor  series), 

-  (1  -  (rA/rMR)2JtMR/(2»)  .  ^  (23) 

so  that  P  in  equation  (14)  is  nearly  linear  with  N. 

P  «  1  -  0.0627N  (24) 
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V . 


For  comparison  purposes,  we  shall  consider  uniform  spacing  where  the  increment 
size  is  independent  of  n.  In  this  case, 

Ar(N)  "  <rMR  "  rA>/N  (25) 

so  that  the  proportion  (14)  does  not  go  negative  until  N  >  26.  However,  it  still 
goes  negative  with  N,  and  this  fact  raises  questions  about  a  discrete  formulation 
of  continuous  phenomena  that  deteriorates  as  the  number  of  elements  becomes  heuristi- 
callv  sufficient  to  replicate  the  continuum.  It  should  be  noted  that  when  N  >  26, 
or  with  the  equal-annuli  algorithm  when  N  >  16,  the  (N  -  l)8^  segment  and  possibly 
other  segments  should  actually  enter  into  the  computation,  and  the  Nth  segment 
should  be  virtually  Ignored. 

The  blade  terminus  for  lifting  purposes  may  be  computed  from  the  summation  of 
the  cuff,  N  -  1  elements,  and  the  scaled  outboard  element: 


N-l 

rT  *  rA  +  £  Arn  +  P  ArN  =  rMRBMR  (26) 

n*  1 


This  formulation  is  independent  of  the  spacing  algorithm.  The  radius  r^  consti¬ 
tutes  the  terminus  of  integration  for  terms  involving  the  lift  coefficient  Cly»  and 
the  suggested  tip-loss  technique  follows: 


Lift 


f  f(r)dr 

*r 

/ 


*r^ 

'rA 

*rN<-) 

"rA 


j  fun 

f (r)dr  +  P  J  f(rN)dr 

rN(-) 


(suggested) 


(extant) 


(27) 


Drag 


JrMR 


g(r)dr 


/ 


rA 

rT 


g(r)dr  + 


|*rMR 


g(rN,rN-l)dr 


(extant) 


(suggested) 


i 


(28) 


With  the  suggested  formulation  the  interval  of  numerical  integration  becomes  indepen¬ 
dent  of  the  size  of  the  Nth  segment,  which  begins  at  rij(_)  and  ends  at  rjj(+)»  so 
that  previous  objections  are  invalidated. 


The  concept  of  a  pseudoblade  is  such  that  it  encompasses  just  the  constant¬ 
lifting-surface  interval  (r^,  r-j).  This  concept  will  be  seen  to  benefit  not  only 
this  analysis  and  the  establishment  of  a  basis,  but  will  also  benefit  the  extant 
simulation  model.  The  technique  uses  linear  extrapolation  on  the  drag  coefficient 
beyond  r^  rather  than  interpolation  (assuming  a  constant  segment  value)  on  the  lift 
coefficient,  and  the  number  N  may  be  as  large  or  small  as  is  desired  without  math¬ 
ematical  inconsistency. 
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Advantages  are  gained  with  the  suggested  technique.  Figure  5  shows  that  lift- 
force  gradients  are  generally  large  near  the  tip,  and  the  suggested  technique  does 
not  approximate  in  this  region.  The  linear  extrapolation  of  drag  over  the  small 
interval  (rj,  r^®)  will  be  shown  to  be  equivalent  to  superimposing  a  perturbation; 
the  magnitude  of  the  lift  profile  in  the  region  of  the  outboard  segment  is  gen¬ 
erally  an  order  of  magnitude  larger  than  that  of  the  drag  profile.  Hence,  lift  space 
is  inferior  to  drag  space  as  an  approximation  medium.  A  final  argument  for  the 
rejection  of  the  extant  technique  concerns  the  distortion  of  continuity,  which  is 
reserved  until  the  end  of  this  section. 


In  order  to  appreciate  the  tip-loss  phenomenon  in  the  context  of  a  pseudoblade, 
we  shall  consider  the  aerodynamic-force  differentials  in  terms  of  the  dimensional 
velocities  in  the  nth  segment  (Houck  et  al.,  1977), 


dFPn  *  \  pcvYn t  (cLYnvTn) /  I cos  Ynl  +  CDYnVPnl^n^ 
dFTn 1 2  3  \  pcvYn<cDYnvTn  '  CLYnVPnk08  Yn|)drn  J 


(29) 


dFRn  *  <vRn/vTn>dFTn 


which  must  be  integrated  over  appropriate  radial  limits  to  obtain  the  blade  aerody¬ 
namic  forces  and  moments.  For  terms  involving  CLy  these  limits  are  the  pseudoblade 
interval  (rA,  rf),  whereas  for  terms  involving  Cpy  these  limits  are  the  total  effec¬ 
tive  blade  interval  (rA,  r^R) ,  which  is  identical  to  the  combination  (rA,  ry)  and 
(rT>  rj®).  This  latter  interval,  defined  as  the  drag  extrapolation  distance,  is 
given  by 


rMR  “  rT  *  (1  ~  ®MR^rMR  =  0,93  ft 


(30) 


which  is  less  than  4%  of  the  total  effective  blade.  This  region  is  always  beyond  the 
outboard  evaluation  point  even  if  N  approaches  infinity.  During  this  small  inter¬ 
val  (rp,  rj®),  the  lift  coefficient  Cly  is  zero  and  the  drag  coefficient  CDy 
along  with  its  appropriate  kinematic  functionality  mey  be  extrapolated  in  both  force 
and  moment  space.  This  procedure  is  outlined  below  in  terms  of  superimposed 
increments : 


AFP (TL) 
AFT(TL) 
AFR(TL) 


1  frMR 

2  PC  J  vyVpCpy  dr 

1  f  rMR 

2  Pc  J  vYvTCDy  dr 

rT 

1  f  rMR 

2  pc  J  vyvRCDY  dr 

rT 


(31) 


Also, 


the  moment  differentials  about  the  hinge  are  defined 


(Houck  et  al.. 


1977)  as 
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1 


dMFA  =  <r  ~  e)dFp  | 

-  (r  ~  e)dFT  j 

so  that  the  extrapolated  moment  Increments  are  given  by 


1 

AMFA(TL)  =  2  PC  J  vYvP^r  "  e)cr 


1  fr> 
AMLA(TL)  *  2  pC  J 

w  v _ 


vyvT(r  -  e)CDY  dr 


By  assuming  linearity  between  the  available  combinations  of  terns  VYnvpnCoYn> 
vYnvTn^DYn’  and  vYnvRn^DYn  at  two  outboard  evaluation  radii  rjj_!  and  rjj,  we 
can  define  four  constant  coefficients: 


rMR 

rT 

rN  “ 

rN- 1 

rMR 

-  rT 

rN  - 

rN-l 

rMR 

-  rT 

rN  “ 

rN-l 

rMR 

"  rT 

rN  - 

rN-l 

<rMR  +  r 


t>  _  rN-l] 


(rMR  +  rMRrT  +  rT>/3]  "  eAF(N-l) 


rwR  —  r«p  r  T 

aM(N)  =  rN  -  rN_:  [(rMR  +  rMRrT  +  rf)/3  -  y  +  rT)J  -  eAF(N) 


so  that  by  using  the  coefficients 

%-l  =  2  pcvY(N-l)CDY(N-l)AF(N-l) 

%  =  2  pcvY (N)CDY (N) AF (N) 

► 

LN- 1  =  I  pcvY(N-l)CDY(N-l)AM(N-l) 

LN  =  2  pCVY(N)CDY(N)AM(N) 

the  indicated  integrations  of  equations  (31)  and  (33)  are  approximated  by: 


I 
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AFP(TL)  =  KN-1vP(N-1) 
aft(tl)  =  KN-1VT(N-1) 
afr(tl)  ”  ^-^(N-l) 
amfa(tl)  =  LN-  1vP  (N-l  > 
amla(tl)  =  ln-ivt(n-i) 


+  %VP(N)' 
+  KNvT(N) 
+  %VR(N) 
+  LNvP(N) 
+  **NVT  (N)  , 


(36) 


The  total  aerodynamic  forces  and  moments  are  thus  generated  by  integrals  of 
equations  (29)  over  the  constant  interval  (rA,  r^)  plus  the  superimposed  tip-loss 
compensation  terms,  that  is. 


FP 


Ff 


F 


R 


M. 


‘FA 


Mi 


LA 


dFp  +  AFp  (XL) 
dFT  +  AFT(TL) 

dFR  +  AFR(TL) 

(r  -  e)dFp  +  AMFA(TL) 

(r  -  e)dFT  + 


(37) 


It  should  be  noticed  that  these  "superimposed  tip-loss  terms”  are  not  related 
to  the  usual  definition  for  tip-loss  effects,  and  their  magnitudes  are  considerably 
less  than  usual  values  because  they  utilize  drag  rather  than  lift  space  as  an  approx¬ 
imation  medium.  Also,  the  linear-extrapolation  process  to  obtain  these  terms  does 
not  influence  the  computational  workload. 


In  the  form  of  equations  (37)  any  continuity  in  the  force  profiles  of  equa¬ 
tions  (1)  may  be  exploited  by  both  the  spacing  and  summation  algorithms.  Since 
these  force  profiles  are  exactly  the  content  of  figures  5  and  6,  reference  to  these 
curves  shows  that  the  physics  of  a  rotor  blade  contribute  enough  discrete  modeling 
problems  without  the  arbitrary  introduction  of  the  tip-loss  discontinuity  into 
the  outboard  evaluation  value.  When  equation  (13)  is  applied  to  these  curves  as 
in  the  extant  model,  the  outboard  evaluation  value  on  both  the  lift  and  drag  pro¬ 
files  usually  is  distorted  to  a  position  that  bears  no  recognizable  relationship  to 
the  other  points,  easily  being  less  than  half  the  magnitude  of  its  adjacent  (N  -  1) 
value,  thereby  creating  a  computer-science  problem  where  none  should  exist.  This 
discontinuous  phenomenon  is  observed  in  all  flight  regimes  and  with  various  values 
of  N.  When  N  ■  5,  for  instance,  Ar^  -  3.12  ft  with  equal-annuli  spacing  and  the 
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Nth  segment  lift  evaluation  is  distorted  to  approximately  70%  of  its  "true"  or  con¬ 
tinuous  value.  Drag  distortion  depends  upon  the  relative  contribution  of  the  terms 
comprising  FT. 

In  the  general  form  of  equations  (37)  the  tip-loss  phenomenon  is  superimposed 
after  the  exacting  tasks  of  the  spacing  and  summation  algorithms  have  been  accom¬ 
plished.  The  suppression  of  distortion  in  the  continuous-derivative  information  is 
extremely  important  when  the  integration  of  equations  (37)  is  approximated  by  a  sum¬ 
mation  process.  Distortion  of  derivatives  should  be  avoided  when  using  discrete 
integration  schemes;  if  any  correlation  is  assumed  between  points,  which  is  certainly 
an  objective  of  discrete  modeling  of  continuous  phenomena  and  is  the  primary  basis 
for  more  sophisticated  integration  schemes,  this  distortion  becomes  computationally 
destructive. 

For  the  20-segments-per-blade  case  the  summed  aerodynamic  lift,  drag,  and  radial 
forces  on  each  blade  ( excluding  the  superimposed  tip-loss  terms  developed  here)  are 
presented  in  figures  9(a)  through  9(c),  and  the  flapping  and  lagging  moments  are 
given  in  figures  9(d)  and  9(e).  These  forces,  shown  for  every  10  knots  of  vehicle 
velocity,  represent  an  accurate  summation  process  operating  on  the  specific  content 
of  figures  5  and  6. 

When  the  tip-loss  ef  feats  are  superimposed,  the  curve  shapes  of  figures  9(a) 
through  9(e)  are  modified  to  become  those  of  figures  10(a)  through  10(e).  It  is 
seen  that  lift-related  terms  vary  negligibly  when  this  tip-loss  technique  is  used 
because  of  the  unique  summation  limits.  Drag  variations  are  small  because  the 
extrapolation  distance  (eq.  (30))  is  small. 

The  quantities  involving  the  multiplier  vp  (perpendicular  velocity)  are  negli¬ 
gible,  as  can  be  seen  by  comparing  the  curves  of  figures  9(a)  and  9(d)  with  those  of 
figures  10(a)  and  10(d),  respectively.  However,  the  quantities  involving  the  multi¬ 
plier  vp  (tangential  velocity)  can  become  significant,  as  can  be  seen  for  the 
advancing  blades  by  comparing  the  drag  curves  of  figures  9(b)  and  10(b).  These 
results  have  been  numerically  investigated  and  the  quantities  AFp(TL)  and  AMjw(TL) 
may  be  arbitrarily  set  equal  to  zero. 

Since  drag  as  a  numerical  quantity  is  usually  smaller  than  lift  and  the  extrap¬ 
olation  distance  is  a  small  and  constant  portion  (0.03)  of  the  total  radius,  Cly 
integrated  to  a  good  approximation  in  this  formulation  from  the  cuff  to  the  aerody¬ 
namic  terminus  for  lifting  purposes  (r-p)  with  a  significant  gain  in  accuracy  (as 
compared  with  the  extant  technique) ,  and  is  integrated  to  the  blade  tip  (r^) 

with  a  slight  degradation  of  accuracy. 

Equations  (34)  through  (37)  constitute  the  superimposed  tip-loss  correction 
terms  that  are  assumed  throughout  the  remainder  of  this  paper. 


V.  SIMPLIFIED  SPACING  AND  SUMMATION 


Three  spacing/summation  algorithms  are  here  examined  in  order  to  reveal  their 
relative-error  convergence  properties  and  to  investigate  their  efficiency  in  the 
rotorcraft  environment.  A  simplified  model  for  force  profiles  is  used  to  demonstrate 
that  the  moment-summation  algorithm  of  the  extant  model,  for  example,  has  slow  con¬ 
vergence  properties.  This  algorithm  is  later  expanded  to  the  full  model  and  account¬ 
able  error  sources  are  identified  in  both  forces  and  moments. 
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In  the  extant  model  the  increment  distribution  is  asymmetrical  about  the  indi¬ 
vidual  radial  evaluation  points  because  of  equal-annuli  spacing,  and  the  summation 
processes  are  given  by  the  quadrature  technique 

N  N 

F  =  rMR  Z  f<yn)Ayn  =  T,  f<rn>Arn 
n=l  n=i 

(38 

N  N 

M=  4  Z  ynf(yn>Ay n  =  E  <rn  "  e>f(rn>Arn 

n=  1  n= 1 

where  yn  is  designated  as  the  center  of  lift-  of  the  segment.  Although  the  limits 
are  proper  in  these  expressions  as  N  approaches  infinity,  it  will  be  shown  that 
they  are  first  approximations  with  relatively  poor  convergence  rates  when  asymmetri¬ 
cal  spacing  is  used. 


In  order  to  demonstrate  convergence  features  of  summation  processes,  a  simpli¬ 
fied  blade-force  profile  is  used;  it  is  developed  as  follows: 


Let  us  consider  a  constant,  unitary  force  where  y  is  the  nondimensional  radial 
variable,  and  neglect  spar  and  offset,  which  complicate  (and  modify  somewhat)  this 
analysis.  Under  these  conditions  the  force-summation  problem  is  to  find  an  efficient, 
discrete  algorithm  that  closely  approximates  the  closed-form  solutions  for  the  force 
and  moment: 


F  « 


M  = 


dy  *  1 

.  1 
y  dy  =  2 


If  this  simplified  model  is  generalized  to  a  force  of  order  k  in  the  indepen¬ 
dent  variable,  the  continuous  solutions  are  given  by 


and  the  direct  summations  (38)  operating  on  this  simplified  force  profile  produce 
relative  force  and  moment  errors  given  by 
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£ 


9  N 

E(F(yk)]  =  (k  +  1)  £ 
n=l 

N 

E[M(yk)]  =  (k  +  2)  £ 
n=l 


k  A 

yn  Ayn 


k+1 

yn 


Ayn  -  1 


1 


(39) 


Since  E[F(yk)]  =  E[M(yk-*)],  the  force  errors  are  subsets  of  the  moment  errors  in 
this  simplified  model  (using  the  summation  piocess  given),  and  they  need  not  be 
included  in  this  discussion.  For  this  reason  the  " relative  moment  error  due  to  a 
fora ^  profile  of  order  k "  will  be  used  to  demonstrate  the  s low-convergence  feature 
of  tft^  extant  moment  summation  algorithm  for  equal-annuli  spacing.  First,  however, 
.a  very  simple  spacing  scheme  is  examined  for  comparison  purposes. 


Uniform-Segment  Algorithm 

We  shall  consider  the  summation  of  moments  for  uniform  spacing  of  segments, 
where  the  force-profile  evaluations  are  made  at  the  midpoints  of  the  segments.  The 
simplified,  nondimensional  radial  parameters  are  given  by 


Ay (n)  =  1/N 

yn=  (n--|)dy(n)  =  (n  -  |)/N 

This  algorithm  will  be  useful  in  isolating  the  error  source  in  the  extant  model. 

The  evaluation  radii  for  various  numbers  of  segments  are  presented  in  table  2. 
(Coincidence  of  the  observation  radii  for  all  lift  and  drag  forces  and  moments  is  an 
assumed  feature  of  a  real  time  model.) 

For  this  algorithm,  as  the  number  of  evaluation  points  is  decreased,  the  end¬ 
point  evaluation  radius  moves  away  from  the  tip  and  the  density  of  evaluations  near 
the  tip  drops  off  in  equal  proportion  to  the  density  elsewhere. 

For  regular  polynomial  forms,  let  us  consider  the  relative  moment  error  from 
equations  (39): 


EIM<^»Un  ■  £  (n  -  IP  -  1  (‘1 

N  n=l 

For  a  force  characteristic  that  is  either  constant,  linear,  or  quadratic,  these 
sums  are 

E[M(y°)]Un  *  0 
E[M(yl)]Un  =  _1/4nZ 
E[M(y2)]Un  =  -1/2N2 

and  they  are  presented  in  table  3  for  up  to  10  evaluation  radii. 
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For  polynomial  forms  this  table  shows  the  general  underprediction  of  moments 
when  the  uniform-segment  algorithm  is  used.  A  priori  information  is  necessary  to 
use  the  simplied  correction  factors.  Factors  such  as  these,  which  imply  an  assump¬ 
tion  about  the  behavior  of  the  force  profiles  between  observation  radii  yn,  are 
later  described  in  terms  of  modifying  the  summation  algorithm  (38).  When  observation 
spacing  is  nor. uniform,  the  effect  under  summation  of  products  such  as  f(yn)Ayn  is 
equivalent  to  a  weighting  algorithm  capable  of  accommodating  an  assumed  curve  profile. 
When  physical  information  is  available  that  indicates  a  tendency  under  certain  condi¬ 
tions  for  some  curves  (such  as  an  advancing-blade  drag  profile)  to  assume  a  particu¬ 
lar  polynomial  form,  this  motivates  the  selection  of  a  nonuniform-spacing  algorithm, 
but  does  not  release  the  spacing/summation  algorithm  from  the  responsibility  for 
accommodation  of  either  simple  or  arbitrary  curve  shapes,  which  may  in  fact  be  the 
dominant  characteristics. 


Equal-Annuli  Algorithm 

In  contrast  to  the  uniform-segment  algorithm,  the  simplified  nondimens ional 
radial  parameters  for  the  equal-annuli  algorithm  (Houck  et  al.,  1977),  which  are 
developed  in  section  VII,  are  given  by 


The  asymmetrical  property  of  these  relationships  is  due  lo  the  fact  that  yn,  the 
evaluation  radius,  is  not  the  midpoint  of  Ayn,  which  is  variable.  For  various 
numbers  of  segments  the  evaluation  radii  are  presented  in  table  4,  which  may  be  com¬ 
pared  with  table  2. 

Unlike  the  uniform-segment  algorithm  of  the  last  section,  this  spacing  algo¬ 
rithm  shifts  the  observation  radii  toward  the  blade  tip.  However,  it  is  also  true 
that  the  endpoint  retreats  from  the  tip  with  a  decreasing  number  of  segments,  which 
is  pertinent  to  the  tip-loss  discussion  of  section  IV. 

For  regular-polynomial  forms  using  the  equal-annuli  algorithm  and  the  direct 
summation  (38),  the  relative  moment  error  for  a  force  profile  of  order  k  is: 


E[M(yk)l 


k  +  2 
( k+2 ) / 2 


(--!) 


(k+i)/z 


(t^n  -  /n  -  1)  - 


By  using  the  Euler-McLaurin  Sum-Formula  (Selby,  1970),  that  is, 

V'1y”=^/'2  1  ,  _i_  1  .  1  11  , 

*-> ,  \3  N  "  2  24N  1920N3  +  9216NS  *  163840N7 

n=  l 


♦•he  series  expansion 
-  0.207886225  ... 


(44) 


the  relative  moment  error  for  a  linear-force  profile  may  be  evaluated, 
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E[M(yl)]Eq  =  W*  E  (n  ~  -  1/11  "  X>  ~  1 

n=  1 


=  0.623658675  ...  -  +  1 


8N2  640N3  "  3072N5 


rT  +  ••• 


(45) 


By  performing  an  absolute  ratio  of  the  relative  moment  errors  the  performance 
of  this  algorithm  relative  to  the  uniform-segment  algorithm  with  k  =  1  is  then 
found 


EfM(y1)]Eq 

ElM(yl)]Un 


2.494634  vfi  -  j  + 


1  1 
160N  “  768N3 


+  .  .  • 


(46) 


which  reveals  that,  for  a  linear  force  and  five  segments,  the  relative  convergence 
rate  of  the  moment  sum  using  the  equal-annuli  algorithm  is  over  five  times  slower 
than  the  moment  sum  using  the  uniform-segment  algorithm.  For  comparison  with  the 
uniform-segment  algorithm  of  table  3,  the  relative  moment  errors  are  presented  in 
table  5  for  the  equal-annuli  algorithm.  Moment  errors  produced  by  using  the  equal- 
annuli  algorithm  thus  compare  poorly  with  those  produced  by  using  the  uniform-segment 
algorithm  of  table  3  when  the  force  profile  is  relatively  simple.  The  large  error 
present  when  the  force  profile  is  very  simple  (e.g.,  constant)  is  disturbing  because 
it  indicates  that  the  extant  simulation  model  produces  avoidable  errors  at  low 
vehicle  velocity.  This  problem  is  isolated  in  section  VII  to  an  inconsistency 
between  the  selection  of  the  interval  Ayn  and  its  observation  radius  yn,  and  this 
problem  is  compounded  when  the  moment  arm  yn  is  used  for  moment  calculations,  as 
above. 


Equal-Annuli  Midpoints 

If  equal-annuli  segment  spacing  is  specified  from  other  consideration...  the  mid¬ 
point  of  the  interval  may  optionally  be  selected  as  the  observation  rad',.*,  i.hat  j  r 


yn  =  (f*n  +  /n  -  1)//n  (47) 

which  preserves  the  spacing  features  of  pure  equal-annuli  spacing  as  shown  in  table  6, 
although  it  gives  slightly  more  weight  to  inboard  segments. 


For  regular  polynomial  forms  using  this  spacing  algorithm  and  the  direct  summa¬ 
tion  algorithm  (38),  the  relative  moment  error  for  a  force  profile  of  order  k  is 
given  from  equations  (39), 


Eimyk)^ 


k  +  2 

2k+lN(k+2)/2 


N 

E 


n=l 


<4T  +  *^7)k 


1 


(48) 


Compared  with  the  pure  equal-annuli  algorithm  shown  in  table  5,  this  algorithm  tends 
to  eliminate  the  large  moment  errors  that  are  observed  for  simple  force  profiles,  as 
shown  in  table  7.  The  constant  column  is  noteworthy.  This  table  indicates,  however, 
that  when  the  force  profile  becomes  quadratic  or  higher  order,  the  gain  in  accuracy 
is  lost.  This  loss  of  accuracy  will  be  shown  in  section  VI  to  be  a  consequence  of 
the  summation  algorithm  rather  than  the  spacing  algorithm. 
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The  use  of  the  simplified  spacing  and  summation  models  for  illustrative  purposes 
has  been  examined,  and  the  following  observations  may  be  made  about  the  summation 
scheme  (38) ,  or  quadrature  technique: 

1.  It  works  well  for  uniform  spacing  and  produces  errors  compatible  with  the 
difference  between  linearity  and  the  actual  order  of  the  force  profile;  that  is, 
table  3  is  quite  elementary. 

2.  It  does  not  work  well  at  all  for  equal-annuli  spacing,  although  some  improve¬ 
ment  is  noticed  if  the  force  profile  is  high  order. 

3.  It  works  better  for  equal-annuli  midpoint  spacing  but  seems  to  deteriorate 
too  rapidly  as  the  curve  order  increases. 

Momentum  or  statistical-sampling  theories  indicate  that  equal-annuli  spacing 
should  deliver  equivalent  or  less  error  than  uniform  spacing,  even  if  the  actual 
force  profiles  are  simple  polynomial  forms,  and  these  approaches  should  all  deliver 
zero  error  when  the  force  profile  is  constant.  The  discrepancy  in  the  extant  model 
is  explored  in  section  VI. 


VI .  SUMMATION  PROCESSES 


The  summation  processes  of  Houck  et  al.  (1977),  where  rn(+)  and  rn(-)  are  the 
upper  and  lower  nth  segment  boundaries,  are  equivalent  in  dimensional  space  to 

N 

F  =  D  f<rnHrn(+)  -  rn(-)J  («> 

n=l 

N 

M  =  £  <rn  ”  e)f(rn)[rn(+)  -  rn(_)]  (50) 

n=  1 

and  have  been  shown  (sec.  V)  to  be  consistent  with  midpoint  spacing  of  the  evaluation 
radii,  that  is. 


rn  =  7  (rn(+)  +  rn(-)l  <51> 

for  elementary  forms  of  the  force  profile  f(r).  However,  they  are  inconsistent  with 
an  asymmetrical-spacing  algorithm  such  as  equal  annuli,  in  which  the  evaluation  radii 
will  be  shown  to  reduce  to 


r 


n 


(52) 


in  that  they  cause  errors  in  both  the  force  and  moment  summations  when  the  force 
profiles  are  elementary  (e.g.,  linear).  This  point  is  important  since  the  very 
nature  of  equal-annuli  spacing  permits  the  acquisition  of  a  limited  amount  of  data 
over  significant  intervals,  that  is,  reduces  the  individual  segment  contributions  to 
what  appears  to  be  elementary  forms  under  summation. 
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When  the  force  profile  is  composed  of  N  perfectly  distributed  (segment 
boundaries  coincident  with  discontinuities)  constant  steps  over  each  (rn(-),  rn(+)), 
the  summation  processes  may  be  examined  at  the  elemental  level.  For  the  force- 
summation  process,  the  force  f(rn)  could  be  evaluated  anywhere  in  the  interval 
(rn(_)  i  ru  <  rn(+))  because  it  is  constant: 

rn  (+) 

f(r)dr  =  f(rn)[rn(+)  -  rn(_)]  (53) 

n(-) 

However,  the  moment-summation  process  (50)  then  requires  rn  to  be  the  midpoint  of 
the  interval. 


/  <r  ~  e)f(r)dr  =  [rn(+)  -  rn(_)]4  [rn(+)  +  rn(_}]  -  e)f(rn)  (54) 

rn (-)  ' 


because  each  element  of  equation  (54)  is  equal  to  each  element  of  equation  (50)  if 
and  only  if  equation  (51)  holds.  This  fact  explains  the  large  errors  observed  in 
column  one  of  table  5. 

When  the  force  profile  is  composed  of  N  perfectly  distributed  (coincident 
boundaries  and  piecewise  discontinuities)  linear  elements,  the  summation  process  may 
again  be  examined  at  the  elemental  level.  In  this  case  the  force  summation  process 
itself  defines  rn  as  the  midpoint.  By  fitting  a  linear  curve  between  the  segment 
extremities, 


rrn(+) 

J  f(r)dr 

rn(-) 


/  n<+>  (ftrn(-)l  +  <ftrn(+)]  “  ftrn(-)l} 
rn(-)  1 


=  1  trn(+)  ~  rn(-)l{f[rn(+)]  +  f[rn(_)]}  (55) 

it  is  shown  that  for  the  force  to  be  linear  its  midpoint  value  must  be  the  average  in 
the  interval, 

f(rn)  =  i  (f[rn(+)]  +  f[rn(_)]}  (56) 

so  that  equation  (51)  is  again  dictated.  Errors  begin  creeping  into  the  midpoint 
algorithm  beyond  this  point  of  complexity,  and  this  very  fact  is  the  motivation  for 
the  selection  of  a  spacing  algorithm  that  tends  to  concentrate  observations  in 
regions  of  high  gradient  activity. 

The  inconsistency  in  the  equal-annuli  algorithm  of  the  extant  model  will  be 
shown  to  be  due  to  the  failure  of  the  summation  processes  (49)  and  (50)  to  accommo¬ 
date  the  spacing  of  equation  (52).  It  should  be  noted  that  equation  (51)  is  con¬ 
sistent  with  these  summation  processes  and  appears  to  be  a  viable  substitute.  For 
this  substitute  spacing  algorithm  the  segment  boundaries  themselves  will  be  shown  in 
section  VII  to  have  an  independent  derivation  based  upon  equal-annuli  arguments. 
However,  objections  still  exist  to  using  equations  (49)  and  (50)  even  if  the  mid¬ 
point  spacing  of  equation  (51)  is  used.  For  instance,  since  the  initial  evaluation 
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radius  rj  is  quite  distant  from  the  cuff  r^  for  both  pure  equal-annuli  and  mid¬ 
point  spacing,  and  since  high  gradient  activity  near  the  tip  is  obvious,  then  equa¬ 
tions  (49)  and  (50)  produce  errors  even  for  piecewise  linear  curves  because  they  fail 
to  extrapolate  exterior  to  rN  (by  using  f(rN_x))  or  interior  to  rj  (by  using 
f(r2))»  (The  interior  slopes  have  a  canceling  action  when  piecewise  linear  curves 
are  used.)  In  addition,  again  assuming  this  linear  form,  the  moment  summation  (50) 
produces  rather  large  errors  because  the  midpoint  radial  arm  is  a  gross  approximation 
to  the  actual  integrated  functionality;  that  is, 

frn(+) 

J  (r  -  e)f (r)dr  *  (rn  -  e)f (rn) [rn(+)  -  rn(_))  (57) 

rn(-) 

An  integration  scheme  that  incorporates  both  the  extrapolation  feature  and  piecewise 
linear  functionality  is  developed  in  the  next  section. 


Integration  of  linear  segments 


In  this  section  a  consistent  summation  algorithm  is  developed.  The  force  pro¬ 
files  are  assumed  to  be  linear  between  observation  radii,  which  are  arbitrary. 


fn(r)  “  fn  +  (rn+1  -  (fn+l  ~  fn>  (1  *  n 


<  N) 


(58) 


and  (if  necessary)  linearly  extrapolated  beyond  the  evaluation  extrema.  With  these 
assumptions,  the  integrations  for  the  forces  and  moments  over  the  entire  pseudoblade 
interval  (r^,  r-j)  may  be  given  by 


F  = 


/ 


r2 


f  1  (r)dr  + 


rA 

r2 


-rT 

"~2  rr  iH-l 

J 

fN-1  (r)dr  +  £  I  f, 

rN-l 

n-2  Jtn 

t  rT 

I^I,2  ^rn+l 

J 

rfN-i(r>dr  +  E  J 

rN- 1 

n=2  rn 

fn(r)dr 


(59) 


(60) 


(excluding  the  superimposed  tip-loss  effects).  The  discrete  solutions  are  given  by 
linear  combinations  of  the  force  observations  fn. 


F  -  E  Vr 


n»l 


«  ■  E  bn£n 
n=  1 


(61) 


when  the  precomputed  coefficients  are  as  follows  for  N  >  3: 
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al  =  2  (r2  ■  rA^2/(r2  -  rl) 

a2  =  \  f r3  ~  rl  ~  (rl  -  rA)2/(r2  ~  rl)l 

aN-l  ”  J  ^rN  '  rN-2  ~  ^rT  ~  rN)2/(rN  “  rN-l^ 

aN  =  J  ^rT  "  rN-l)2/(rN  “  rN-l) 

an  ~  ~2  ^rn+i  ~  rn-i^  (2  <  n  <  N  -  1) 

and 

bl  =  \  (r2  -  rA) 2  (r2  +  2tp)K*2  “  rl>  “  eai 

b2  =  -|  [ (r3  ~  r2> (r 3  +  2r2)(r2  -  rx)  +  2 (r|  -  r|) 

-  3r  i  (r2  -  r^)]/(r2  -  rx)  -  ea2 

bN-l  “  5"  KrN-l  -  rN-2^rN-2  +  ^rN-l^rN  ”  rN-p  +  ^rN^rT  “  rN-l^  " 

“  2<rT  “  rN-l)]/<rN  "  rN-l>  "  eaN-l 

bN  =  6'  (rT  “  rN-l^2(rN-l  +  2rT^ /  (rN  ~  rN-l^  _  eaN 

bn  =  \  <rn+l  "  rn-l)<rn+l  +  rn  +  rn-l>  “  ean  (2  <  n  <  N  -  1) 

- 


(62) 


(63) 


Hence,  by  merely  assuming  a  linear  transition  of  the  observed  force  profiles  from 
point  to  point,  a  consistent  set  of  coefficients  are  made  available  for  general  sum¬ 
mation  application  by  using  equations  (61);  they  accommodate  any  monotonic  spacing 
algorithm  rA  -  rk  <  r^-n  ^  rj  (1  £  k  <  N),  and  they  extrapolate  internal  to  the 
inboard  and  external  to  the  outboard  observation  radii.  In  addition,  for  moment 
summation  a  linear  force  profile  is  error  free. 


These  general  summation  coefficients  are  independent  of  the  spacing  algorithm, 
and  they  will  be  applied  to  the  basis  in  the  next  section.  Later,  they  will  also  be 
applied  to  the  "equal-annuli  midpoint  algorithm"  with  accuracy  improvement  over  the 
extant  model. 


When  N  *  3,  the ‘summations  in  equation  (59)  vanish,  so  that: 
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a2  =  j  [r3  -  rx  -  (r3  -  rA)2/(r2  -  -  (rT  -  r3)2/(r3  -  r2)] 

b2  =  I  (r^  1  r i)  [~3rl(r2  +  rA)  +  2 (r|  +  r2rA  +  rA) ] 

+  J  —  r^V3r3(rT  +  r2>  “  2(rT  +  rTr2  +  r2)]  -  ea2 


Algorithm  for  the  Basis 

The  basis  always  utilizes  20  equally  spaced  observation  radii.  The  initial 
evaluation  radius  is  at  the  cuff  rj  =  rA,  and  the  final  evaluation  radius  is  at  the 
tip  for  lifting  purposes  rN  =  r-j.  This  selection  obviates  extrapolation  except  for 
the  tip-loss  correction,  and  it  assures  that  the  best  available  parameters  are  used 
for  the  tip-loss  extrapolation  process.  The  basis  also  uses  the  piecewise  linear 
summation  process  of  the  previous  section.  The  summation  coefficients  reduce  to 

[\  Ar  (n  =  1) 

(1  <  n  <  N)  (64) 

(n  =  N) 


’  i 

A-  Ar(rA  +  Ar/ 3  -  e) 
bn  =  '  Ar<rn  -  e> 

4  Ar(rT  -  Ar/3  -  e) 

\ 

where 

Ar  =  (rT  -  rA)/(N  -  1) \ 

}  (66) 

rn  =  rA  +  (n  -  l)Ar  f 

The  reduction  of  the  coefficients  to  these  simple  forms,  which  appear  similar  to 
those  of  the  extant  model,  is  due  to  the  selection  of  r3  to  correspond  to  the  cuff, 
rfj  to  correspond  to  the  tip  for  lifting  purposes,  and  uniform  segments.  When  a 
nonuniform-spacing  algorithm  is  used  this  simple  form  does  not  occur. 

Because  of  the  occasional  high  rate  of  curvature  change  in  the  outboard  portion 
of  the  drag  profiles,  if  it  is  required  that  the  number  of  segments  be  minimal  ('5), 
this  particular  subset  of  the  linear-integration  algorithm  is  not  recommended.  How¬ 
ever,  as  shown  in  figures  5  and  6,  since  the  plotting  software  has  conveniently 
drawn  straight  lines  through  each  of  the  inclusive  observation  points,  it  is  easily 
seen  that  by  using  20  observations  (standard  for  the  basis)  this  linear-integration 
algorithm  produces  negligible  differences  from  the  continuum. 

Figure  10  has  been  created  by  using  this  20-segment  algorithm  and  including 
the  superimposed  tip-loss  terms  as  given  in  section  IV.  When  resolutions  of  the 


(n  =  1) 

(1  <  n  <  N)  (65) 

(n  =  N) 
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aerodynamic  torces  and  moments  are  performed  and  inertial  terms  are  included,  the 
rotor  system  delivers  body-axes  forces  and  moments  as  given  in  figure  11.  Fig¬ 
ures  10  and  11  thus  constitute  the  standard  for  comparison  that  may  be  used  to  estab¬ 
lish  the  efficiency  of  alternate  models. 

Integration  of  Quadratics 

It  is  maintained  that  if  the  spacing  algorithm  accounts  for  curvature  in  the 
force  profiles  via  the  use  of  unequal  observation  "windows,"  then  the  summation 
algorithm  may  be  proportionally  simplified.  In  particular,  for  the  RSRA  data  the 
linear-segment-integration  technique  discussed  previously  is  sufficient  provided  only 
that  the  observation  radii  are  appropriately  placed  within  each  window.  In  order  to 
demonstrate  this  fact  a  higher-order  integration  process  was  also  examined. 

As  in  the  linear-segment- integration  case,  a  set  of  coefficients  may  be  devel¬ 
oped  that  give  a  consistent  summation  process  when  the  behavior  between  points  is 
assumed  to  be  quadratic.  Of  course,  at  least  three  observations  are  required.  If 
the  observation  radii  are  assumed  to  be  rn,  rn+i,  and  rn+2,  the  pertinent  integrals 
of  an  assumed  quadratic  function  f(r)  over  the  interval  (rx,  ry)  are  given  by  com¬ 
binations  of  the  observed  fn,  where  k  =  0  or  k  =  1: 


_k+l  _k+l\ 
rv  rx  > 


rkf(r)dr  ■  (rn+2  -  rnMrn+1  -  rn)  [k  +  1  rn+2rn+i  (r£+1  -  r£+1) 

-  in ~2  (rn+2  +  rn+l><ry+2  "  rx+2>  +  FT3  (ry+3  “  rx+3)] 

_ !n±i _  f  1  r  r  /  k+1  _  k+l\ 

<'„«  -  '„«><'„«  -  ■•»)  U  +  1  "+2”ly  r*  ’ 

-  iTTT  (rn+2  +  rn)(4+2  -  rk+2)  +  ^  (r£+3  _  rk+3)j 

+ _ _ f  1  r  r  frk+i  _  rk+i\ 

(rn+2  -  -n)(rn+2  -  rn+l>  U  +  1  rn+lrn(ry  rx  > 

-  iTTT  (rn+i  +  rn><ry+2  "  rx+2)  +  kT3  (r£+3  ■  rx+3)]  (67) 

The  range  of  integration  or  applicability  of  equation  (67)  must  be  decided,  and 
this  process  becomes  somewhat  complicated  when  N  is  general.  However,  for  n  =  5, 
a  fairly  natural  selection  seems  to  be  to  pivot  two  quadratics  about  the  observation 
at  r3.  In  this  case  the  range  of  (rx,  ry)  would  be  (r^,  r3)  for  the  first  quadratic 
and  (r3,  rj)  for  the  second.  In  this  case  the  summation  coefficients  for  use  in 
equation  (61)  are: 
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ai  =  £r3r2(r3  -  rA)  -  j  (r3  +  r2)(r3  -  r£)  +  (r3  -  r|)/3j/[(r3  -  ri)(r2  -  r^] 

a2  =  -|^r3r1(r3  -  rA)  -  (r3  +  rx)(r3  -  r£)  +  (r3  -  rl)/3j/[(r3  -  r2)(r2  -  rj)] 

a3  =  [r2rl(r3  -  *A)  "  \  +  rl><r3  “  rA>  +  (r3  ~  ri)/3]/[(r3  -  rj)(r3  -  r2)] 

+  £r5r4(rT  -  r3)  -  j  (r5  +  r4)(r£  -  r3)  +  (r|  -  r3)/3  /[(r5  -  r3)(r4  -  r3)] 

a4  *  _ £r5r3 (rT  "  r3>  “  Y  (^5  +  r3)(rj  -  r3)  +  (r|  -  r3)/3j/[(r5  -  r4)  (r4  -  r3)] 

a5  *  [r4r3<rT  -  r3>  -  \  (r4  +  r3) (rT  -  r 3 )  +  (r|  -  r3)/3j/[(r5  -  r3)(r5  -  r4)] 

(6 

and 

bl  =  [f  r3r2<r3  ~  rA>  “  3  <r3  +  r2)  (r3  ~  r|)  +  i  (r3  -  r£)]/[(r3  -  rjXr;,  -  r^] 

-  eax 

b2  -  r3r l  (r3  “  r\)  -  i  (r 3  +  ri)(r3  -  r^)  +  j  (r^  -  4>]/[(r3  -  r2)  (r2  -  r^] 

-  ea2 

b3  =  [f  r2r i (r3  -  r^)  -  -j  (r2  +  ri)(r|  -  r,|)  +  ^  (r^  -  r£)]/[(r3  -  ri)(r3  -  r2) ] 

+  r5r4(rT  -  r3)  -  <r5  +  r4)(r|  -  r3)  +  (r£  -  r3)  /[(r5  -  r3)  (r4  -  r3)] 

-  ea3 

b4  =  r5r3 <rx  "  r3)  -  -j  (r5  +  r3)(r|  -  r|)  +  j  (r£  -  r^)J/[(r5  -  r4)(r4  -  r3)] 

-  ea4 

b5  ”  r4r  3  (rT  -  r3>  -  3  (r4  +  r3)(r|  -  r3)  +  ^  (r£  -  r3)J/[(r5  -  r3)(r5  -  r4)  ] 

-  ea5  f 


Interestingly  enough,  the  use  of  these  coefficients  produces  negligible  improve¬ 
ment  over  the  use  of  the  linear  coefficients  given  in  equations  (62)  and  (63)  when 
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the  "equal-annuli  midpoints"  algorithm  (developed  in  sec.  VII)  is  used  to  select  the 
observation  radii.  Hence,  higher-order  summation  approximations  such  as  this  are 
superfluous  when  the  spacing  algorithm  itself  accounts  for  the  curvature  of  the 
force  profiles. 


VII.  EQUAL-ANNULI  DERIVATION 


Since  a  more  general  and  accurate  summation  process  has  i.,  k  -led,  we  now 
turn  our  attention  to  the  spacing  algorithm.  We  shall  consider  ..e  nzinn. _ n  advance 
angle  dip  that  occurs  during  some  time  interval  such  that  the  di-  tan  mat  a  point 
on  the  blade  at  the  observation  radius  rn  travels  is 

dSn  =  rn  d<p 

and  the  area  swept  out  from  the  hub  to  this  point  is 

d^n  =  y  rn  dSn  =  j  rn  d^ 

The  equal-annuli  concept  is  that  the  area  swept  out  by  an  adjacent  observation  radius 

dAn+l  =  ~2  rn+l  dsn+l  *  y  rn+l  d^ 

is  used  to  define  the  differential  area 

AA  *  dA^+j  ~  dAn  *  ~2  (fn+1  ~  rn)di(/  (70) 

which  is  independent  of  n.  Hence,  the  difference  in  the  squares  of  radii  is 
constant, 

K  =  4+1  -  4  (71) 

which  is  another  way  of  saying  that  the  density  of  observations  increases  hyperboli- 
cally  with  the  midpoints  between  observations,  that  is, 

I  K 

rn+l  -  rn  =  T - - -  <72> 

2  (rn+l  +  rn> 

By  induction,  equation  (71)  may  be  written  in  terms  of  the  initial  radius  ri- 

r^  =  (n  -  1)K  +  4  (73) 

The  constant  K  is  determined  by  dividing  the  total  area  swept  out  by  the 
effective  blade  (r^  used  here  rather  than  r^)  into  N  equal  areas, 

(r|  -  r*)dip  *  NK  dip  (74) 


27 


so  that 


r 


2 

n 


n  -  1 


N 


rfr  +  r 


2 

1 


(75) 


The  initial  observation  radius  is  obtained  by  noting  that  the  area  from  r^  to  rj 
is  half  the  area  from  rj  to  r2»  which  is  equivalent  to 

r?  -  ri  =  \  (cl  -  rf )  (76) 

and  substituting  this  into  equation  (75)  with  n  =  2  yields 

rl  =  rj  +  (r|  -  rJ)/(2N)  (77) 


so  that  (when  rj^  is  substituted  for  rj)  the  blade  observation  radii  are  given  by 
the  reference  result  (Houck  et  al.,  1977), 


r !  =  [r2  +  (r2  -  r2)/(2N)]1/2 

rn  =  trn-l  +  (rT  -  4)/n]1/2  (1  <  n  £  N) 


(78) 


which  by  induction  may  be  written  for  all 

rn  =  [rA  +  (rT  " 


(79) 


As  a  spacing  algorithm  this  one  appears  satisfactory  because  it  tends  to  appre¬ 
ciate  the  occasional  curvature  change  in  drag  forces  with  radius  or  at  least  recog¬ 
nize  the  linear  radial  factor  required  for  the  computation  of  aerodynamic  moments. 

A  summation  algorithm,  however,  should  recognize  the  convergence  problem  associated 
with  the  approximation  processes. 


mVj  N 

F  =  I  f  (r)dr  =  £  f(rn)an 

rA  n=1 


M  = 


(r  -  e)f(r)dr 


N 


E 


f (rn)bn J 


(80) 


and,  therefore,  the  coefficients  (an,  bn)  should  be  determined  from  assumptions  about 
the  behavior  of  f(r)  between  observation  radii  rather  than  be  averaged  over  seg¬ 
ments  that  are  artifacts  in  the  development  of  the  spacing  algorithm.  This  process 
is  discussed  in  the  next  section. 


Equal-Annuli  and  Quadratures 

In  this  section  the  segment  spaaing  for  the  equal-annuli  algorithm  is  developed. 
It  is  noted  that  the  selection  of  observation  radii  within  these  segments  involves 
other  considerations  when  the  quadrature  technique  is  used  in  the  derivation. 
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Let  us  consider  the  segment  bounded  by  rn(+)  on  its  exterior  and  rn(-)  on  its 
interior.  The  distances  that  these  extremity  positions  travel  during  some  azimuth 
advance  angle  dip  are  given  by 

dSn(+)  =  rn(+)d* 
dSn(->  =  rn(-)d^ 


so  that  the  total  area  swept  out  by  the  segment  is 

AA  =  I  frn(+)  "  rn(-)Jd^ 

When  these  segment  areas  are  independent  of  n  we  have 

K  "  rn (+)  '  rn (-) 

so  that  by  using  the  contiguous  relationships 


rn(+)  -  rn+l (-) 


rn(-)  rn- 1 (+) 


the  indices  may  be  made  compatible: 


K  =  rn (+)  -  rn-l (+) 


K  =  r. 


n+l (-)  “  rn (-) 


The  extremities  of  the  segment  are  then  found  by  induction, 


rn(+) 


=  nK  +  L 


r 1/  x  =  (n  -  1)K  +  L 


where  the  constants  are  determined  from  the  boundaries 


rN(+) 


and  the  solutions  are  given  by: 


rn (+)  =  trA  +  <rT  -  4>n/N31/2 
rn(_)  =  [r£  +  (rx  -  r£) (n  -  1)/N] 


The  segment  length  may  be  expressed  by 
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Arn  =  rn(+)  -  rn(-) 


but  it  has  not  been  relevant  to  this  development  since  only  the  differences  in  the 
squares  of  radii  have  been  used.  Furthermore,  one  may  define  a  radius  such  as  in 
the  extant  algorithm  (eq.  (79)),  but  it  has  not  been  relevant  to  this  development 
either.  A  possible  reason  tor  its  inclusion  in  the  algorithm  for  the  extant  model 
is  as  follows:  If  the  assumption  is  made  that  the  force  between  points  is  quadratic 
in  radius,  specifically  at  the  extremities,  that  is, 

fn (+)  =  Crn(+) 

fn(-)  =  Crn(-) 

then  another  assumption  is  necessary,  which  is  inconsistent  with  the  first  assump¬ 
tion:  that  the  average  of  the  functional  values  at  the  segment  extrema  is  equal  to 

the  observation  value,  that  is, 

fn  =  \  tfn (+)  +  fn(-)l  =  J  C[r£(+)  +  r£(_)]  (88) 

so  that,  by  applying  the  original  assumption  to  the  evaluation  point  itself, 

fn  =  Crn  (89) 


the  simultaneous  solution  of  equations  (88)  and  (89)  produces  the  formulation  of  the 
extant  model,  which  is  equivalent  to  equation  (79): 


(90) 


But  if  f(r)  is  indeed  quadratic  as  in  equation  (89),  then  the  correct  integration 
process  is  given  by 


frn(+)  „  fv  n(+) 

I  f (r)dr  =  C  I  r2 

*M-)  *'rn(-) 


3  C[r^(+)  "  rn(-)1[rn(+)  +  rn (+)rn(~)  +  rn(-)^ 


(91) 

so  that  the  correct  quadrature  method  of  integrating  the  forces  within  the  interval 
is 


f(rnMrn(+)  ”  rn(-)  1  ~  ^rnl-rn(+)  rn(-)^ 

=  3"  C[rn(+)  '  rn(-)l[^(+)  +  rn(+)rn(-)  +  rn(-)l 


(92) 


This  expression  results  in  an  observation-radius  definition  in  terms  of  the  segment 
extremities,  which  is  quite  different  from  the  reference  spacing: 


rn  =  ]/\  lrn(+)  +  rn (+)rn(-)  +  ^(-)1 


(93) 
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Thus,  this  particular  segment-oriented  derivation  of  the  equal -annuli  algorithm 
isolates  the  observed  mathematical  inconsistency.  Equal-annuli  spacing  of  the  ob¬ 
servation  radii  is  not  compatible  with  the  quadrature  technique.  However,  as  given 
in  equation  (79),  equal  annuli  spacing  may  be  used  with  the  summation  coefficients 
of  Section  VI.  Indeed,  this  combination  is  recommended  for  real-time  simulation 
programs,  where  a  minimal  number  of  radial  observations  are  required. 


Equal-Annuli  Midpoints 

The  combination  of  the  quadrature  technique  and  equal-annuli  spacing  tends  to 
improve  results  when  the  force  profile  is  a  high-order  polynomial  in  radius.  How¬ 
ever,  this  feature  is  rarely  of  benefit,2  while  its  destructive  tendencies  are  much 
more  common.  In  order  to  illustrate  this  process,  an  alternate  spacing  algorithm 
is  introduced,  which  is  more  compatible  with  the  quadrature  technique  because  it 
uses  the  midpoints  of  the  artificial  segments,  given  in  equations  (87): 


r 

n 


(rT  “ 


+  r: 


(r2  - 


+  r 


(94) 


For  a  quadratic  force  profile,  this  spacing  algorithm  is  shown  in  appendix  A  to 
deliver  half  as  much  error  as  equal-annuli  spacing,  when  the  quadrature  technique 
is  used. 

The  spacing  algorithm  of  equation  (94)  is  also  of  interest  when  the  suggested 
summation  coefficients  (62)  and  (63)  are  used. 


F  ■  Z 

n=l 


M  =  £  V<rn> 

n=l 


(95) 


In  appendix  B  this  process  is  compared  to  the  original  combination  of  equal-annuli 
spacing  and  quadrature  technique.  The  use  of  summation  coefficients  is  shown  to 
have  a  distinct  advantage  for  arbitrary  curve  shapes. 

In  combination  with  the  tip-loss  technique  of  equation  (37),  the  inconsisten¬ 
cies  in  the  radial  computations  are  eliminated  for  both  real-time  simulation  and 
for  continuum  emulation,  provided  the  equal-annuli  spacing  algorithm  of  equation 
(79)  or  the  alternate  spacing  algorithm  of  equation  (94)  is  used,  along  with  the 
summation  coefficients.  No  real-time  computational  penalty  accrues.  Of  particular 
importance  to  real-time  simulation  (using  a  minimal  number  of  radial  observations) 
are  the  features  of  linearity  between  observation  radii,  extrapolation  internal  to 
the  inboard  observation  radius  and  external  to  the  outboard  observation  radius, 
superimposed  tip  loss  that  does  not  influence  the  continuity  of  the  functionals  to 
be  summed,  a  consistent  solution  for  forces  and  moments,  and  no  extra  real-time 
computational  expense. 

20nly  the  drag  profiles  for  advancing  blades  during  high-speed  flight  benefit 
(e.g.,  blade  2  in  fig.  6 (ee) ) . 
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VIII.  ERROR  COMPARISON 


The  errors  from  both  the  extant  and  alternate  algorithms  are  examined  here  by 
using  5-segment  models  and  comparing  them  with  the  20-segment  basis.  In  order  to 
fully  appreciate  this  comparison  the  "constrained  and  balanced"  operations  are 
necessary. 

In  exercising  the  extant  model,  only  the  angle-of-attack  correction  of  sec¬ 
tion  III  was  implemented. 

In  figure  12  the  total-body-axis  forces  and  moments  for  the  extant  model  are 
presented  in  the  constrained  condition.  Thus,  the  blade  trajectories  have  been  con¬ 
strained  to  the  same  trajectories  as  the  20-segment  basis  case.  The  forces  and 
moments  of  the  basis  of  figure  11  (the  basis  always  uses  20  segments)  are  also 
repeated  on  these  graphs  for  comparison  purposes.  The  differences  are  due  to  the 
decreased  number  of  segments,3  slow  summation  convergence,  the  tip- loss  computation 
technique,  and  equal-annuli  spacing. 

Since  figure  12  has  been  created  in  the  constrained  condition  so  that  inertial 
terms  are  invariant,  differences  are  entirely  caused  by  the  computation  of  aerody¬ 
namic  terms.  These  aerodynamic  errors  are  presented  in  figure  13.  Among  items  that 
should  be  noticed  are  the  force  and  moment  standoff  differences  at  low  vehicle  veloc¬ 
ity,  the  deterioration  (see  blade  5  in  fig.  13(b))  of  reverse-flow  computation,  and 
error  variation  when  the  blade  encounters  drag  divergence  (see  blade  2  in  fig.  13(a)). 
From  the  data  set  this  latter  phenomenon  is  identified  by  variations  in  the  drag 
coefficient  due  to  Mach  number  variation  without  regard  to  angle-of-attack  variation. 

When  the  blade  states  are  allowed  to  interact  with  their  flap/lag  differential 
equations,  the  balanced  condition  is  achieved.  Figure  14  shows  the  extant  and  basis 
rotor-system  total-body-axes  forces  and  moments.  The  new  spatial  orientation  of  the 
rotor  disc  produces  differences,  some  of  which  involve  inertial  terms.  Although  the 
forces  of  figures  13(a)  and  14(a)  do  not  experience  much  variation,  an  interesting 
phenomenon  occurs  in  moment  space,  especially  in  yawing  moment;  it  is  seen  in  com¬ 
paring  the  value  t  R^;  between  figures  12(b)  and  14(b).  In  order  to  explain  this 
phenomenon  figure  15  is  presented,  which  gives  the  rotor-torque  differences.  This 
figure  shows  that  the  torque  difference  under  the  constrained  condition  is  minimal; 
however,  when  the  flapping-  and  lagging-moment  variations  of  figures  13(d)  and  13(e) 
are  allowed  to  interact  with  the  differential  equations,  new  blade  states  are  created 
that  minimally  influence  the  force  values  but  cause  a  rather  large  variation  in  total 
rotor-system  torque.  This  effect  is  shown  in  the  following  section  to  be  attribu¬ 
table  to  one  specific  inertial  term. 

The  alternate  simulation  model  creates  errors  that  are  more  consistent  with  its 
assumptions,  and  usually  much  less  than  those  of  the  extant  model.  Individual-blade 
aerodynamic  force  and  moment  errors  for  this  five-segment  model  are  given  in  fig¬ 
ure  16  for  comparison  with  figure  13.  The  total  system  outputs  in  the  constrained 


3At  low  vehicle  velocity  the  unusually  large  displacement  of  the  initial  evalua¬ 
tion  radius  for  the  extant  algorithm  delivers  an  error  in  perpendicular  force  of 
approximately  100  lb  per  blade.  Increasing  the  number  of  segments  from  5  to  10  only 
reduces  this  error  to  approximately  70  lb  per  blade.  The  extrapolation  feature  of 
the  alternate  model  eliminates  this  error. 
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condition  are  given  in  figure  17  for  comparison  with  the  extant -model  outputs  shown 
in  figure  12.  The  total  system  outputs  in  the  balanced  condition  are  given  in  fig¬ 
ure  18  for  comparison  with  the  extant-model  outputs  shown  in  figure  14.  As  might 
be  expected,  the  total  torque  error  is  much  improved.  This  error,  given  in  figure  19 
for  both  the  balanced  and  constrained  conditions,  may  be  directly  compared  with  the 
error  for  the  extant  model  shown  in  figure  15;  a  further  examination  of  the  torque 
error  is  given  in  the  next  section. 

Figures  1  and  2  have  been  presented  to  demonstrate  thrust  variations  with  seg¬ 
ment  number.  Only  on  these  two  graphs  does  the  "basis"  contain  less  than  20  ele¬ 
ments.  The  improvement  in  thrust  obtained  by  using  the  alternate  algorithm  rather 
than  the  extant  algorithm  is  demonstrated  in  these  figures  for  vehicle  velocities  of 
10  and  250  knots.  Specifically,  for  a  small  number  of  segments  the  alternate  algo¬ 
rithm  is  at  least  twice  as  accurate,  and  for  a  larger  number  N  the  convergence  is 
faster.  Invariancy  with  N  is  displayed  by  the  basis  in  the  region  of  N  =  20. 
Another  feature  that  is  indicated  by  these  graphs  is  that  lift  forces  are  improved 
with  uniform  spacing  (basis) .  Because  of  their  nearly  arbitrary  variations  this 
improvement  with  uniform  spacing  is  most  probably  true  in  general,  but  this  fact  is 
not  true  for  drag  forces.  These  drag  forces  exhibit  distinct  polynomial  forms, 
especially  at  significant  velocities  for  advancing  blades,  and  uniform  spacing  fails 
to  accommodate  this  feature  when  the  number  of  segments  is  small.  Indeed,  high- 
order  polynomial  drag  characteristics  are  the  justification  for  equal-annuli  spacing, 
and  this  feature  has  been  preserved  in  the  development  of  the  alternate  algorithm. 

It  should  be  noted  that  this  alternate  algorithm  is  also  consistent  with  lower-order 
curves,  which  are  dominant  in  figures  5  and  6. 


Error  Propagation 

The  extant  model  is  here  used  to  trace  an  error  from  its  source  in  the  con¬ 
strained  condition  to  the  balanced  condition.  It  will  be  shown  ti  at  the  major  dif¬ 
ferences  in  balanced  and  constrained  conditions  are  simply  due  to  an  inertial  term. 
(This  fact  supports  the  contention  that  a  major  contribution  to  rotorcraft  trimming 
will  result  from  the  "constrained"  concept.)  In  the  previous  section  the  importance 
of  accurate  summation  along  the  blade  was  illustrated  by  a  variation  in  torque 
(fig.  15).  Although  the  aerodynamic  force  outputs  during  the  balanced  condition  vary 
only  slightly  from  the  outputs  during  the  constrained  condition  (as  can  be  seen  by 
comparing  figs.  12(a)  and  14(a)),  the  slight  reorientation  of  the  rotor  disc,  caused 
by  the  inaccurate  moment  summation  process,  produces  torque  errors.  These  errors  are 
due  primarily  to  the  contribution  of  a  single  inertial  term,  as  shown  below. 

The  total  rotcr  system  yawing  moment  is  nearly  identical  to  the  torque  Q. 

Slight  differences  are  attributable  to  the  moment  arm  and  to  the  shaf t-to-body-axis 
transformation,  which  is  almost  an  identity  matrix,  hence  for  B  blades, 

B 

nMR  =  Q  =  "23  (MLDi  cos  Bi  +  eFXAi  +  eFXIi^  (9&) 

i=l 

For  the  purpose  of  a  perturbation  analysis,  we  consider  a  portion  of  the  inertial 
term, 


33 


B  B 

Nwr  -  -e  £  FXIi  +  ...  =  -e  Y  Mb  cos  &±  sin  6±(rs  -  ft)2  +  ... 
i=l  1=1 

B 

=  -eMbft2  Y  6i  +  •••  ~  ~eBMbn2aOL/57.3  +  ...  (97) 

i=l 

where  the  principal  lagging  Fourier  coefficient  is  defined  (in  deg)  as 

a0L  =  ^  E  6i  <98> 

Thus,  if  the  torque  error  is  due  primarily  to  a  reorientation  of  the  lagging  angles, 
it  may  be  reconstructed  by  considering  only  the  contribution 

ANMR  a  _eBMb^2  AaOL/57.3  (99) 

where  rpm  follows  the  computed  schedule  of  table  1. 

The  lagging  coefficient  agL  is  presented  in  figure  20  from  the  basis.  Also 
appearing  on  this  graph  is  the  difference  in  lagging  coefficient  AagL  observed  from 
the  extant  simulation  model  with  balanced  conditions  for  five  segments.  Using  equa¬ 
tion  (99),  with  the  parameters  e  =  1.05,  Mb  =  91,  B  =  5,  and  rpm  as  in  table  1, 
figure  21  is  produced,  showing  very  little  deviation  from  the  torque-error  history 
of  figure  15.  Hence,  inaccuracies  in  the  computation  of  the  lagging  moment  during 
the  constrained  condition  are  seen  in  the  balanced  condition  to  contribute  directly 
to  the  yawing-moment  error,  as  should  be  expected. 

From  this  comparison  an  important  fact  is  revealed:  In  order  to  obtain  accurate 
rotor  system  outputs  the  moment  summations  are  at  least  as  important  as  the  force 
summations.  The  amplification  due  to  the  radial  moment  arm  means  that  any  effort 
spent  in  replicating  the  true  force  profiles  is  effort  well  spent. 

The  lagging-moment  error,  which  directly  influences  the  lag  states  through  the 
forcing  functions  of  the  differential  equations  (Houck  et  al.,  1977),  influences 
agL  via  equation  (98),  and  thus  contributes  to  the  torque  error  via  equation  (99). 

It  has  been  shown  here  that  this  logic  path  is  the  one  that  produces  the  torque  error 
for  the  extant  simulation,  which  manifests  itself  in  the  total-yawing-moment  errors 
given  in  figure  14(b). 


IX.  CONCLUSIONS 


An  a.'gle-of-attack  correction  for  reverse-flow  computation  has  been  developed  by 
using  symmetry  arguments,  and  this  correction  has  been  demonstrated  to  influence 
simulation  results  at  vehicle  speeds  in  excess  of  155  knots. 

The  balanced/constrained  dynamics  concept  has  been  developed  for  rotor  systems 
and  shown  to  be  a  valuable  tool  for  comparing  different  formulations.  It  should  also 
lead  to  the  development  of  a  generalized  rotorcraft  trimming  process  that  promises  to 
be  computationally  efficient. 
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An  alternate  tip-loss  model  has  been  developed  that  simplifies  the  integration 
process  by  using  a  minor  perturbation  that  is  superimposed  upon  integrals  of  contin¬ 
uous  derivatives.  The  pseudoblade  concept  has  been  demonstrated  to  exploit  the  rela¬ 
tive  magnitudes  of  lift  and  drag  forces  for  the  approximation  process. 

The  slow  convergence  properties  of  the  radial  summation  scheme  using  quadratures 
has  been  noted,  and  the  propagation  of  errors  has  been  illustrated.  The  summation 
scheme  has  been  improved  with  a  set  of  coefficients  that  are  applicable  to  various 
spacing  algorithms;  they  are  demonstrated  herein  by  using  the  alternate  algorithm 
that  features  equal-annuli  midpoints.  Equal-annuli  spacing  is  recommended  for  real¬ 
time  use  with  these  summation  coefficients.  Output  improvements  are  shown  to  be 
especially  noteworthy  during  low  vehicle  velocity.  A  general  improvement  in  the 
computation  of  lift-related  quantities  has  been  demonstrated,  although  at  higher 
flight  velocities  some  deterioration  has  been  observed  in  drag-related  quantities. 
Since  lift  quantities  are  usually  an  order  of  magnitude  greater  than  drag  quantities, 
the  total-rotor-system  outputs  are  improved. 

The  advancing-blade  drag  profiles  are  shown  to  occasionally  exhibit  a  high-order 
radial  polynomial  behavior  that  is  best  accommodated  by  a  biased  spacing  algorithm 
such  as  in  the  extant  model,  but  this  advantage  is  deleterious  to  the  convergence  of 
other  quantities,  including  the  retreating-blade  drag  profiles.  Further  research  is 
recommended  in  the  area  of  summation  coefficients  that  recognize  both  the  character¬ 
istic  differences  in  lift  and  drag  and  the  azimuth-dependency  problem. 

A  data  base  of  high  fidelity  has  been  created  by  using  20  uniform  segments. 

This  data  base  has  been  graphically  illustrated  and  exists  as  a  basis  for  additional 
model  development  activities.  This  basis  exhibits  neither  standoff  errors  nor  eval¬ 
uation  weighting  usually  associated  with  real-time  algorithms  that  require  a  minimal 
number  of  evaluations;  it  is  virtually  invariant  with  the  number  of  segments  used. 

For  a  large  number  of  segments,  the  alternate  algorithm  converges  to  the  basis. 

VJhen  tip  loss  is  reformulated  as  developed  herein,  the  extant  algorithm  also  con¬ 
verges  to  the  basis,  but  more  slowly. 
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APPENDIX  A 
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c 


JEq  2 


1 

=  2  ^  [rn(+)  +  rn(-)Hrn(+)  -  rn(_)]  -  Fc 


n=  1 


i  n  _  _ 

=  T  £  f<2n  -  Dbl  +  2b2][/nb1  +  b2  -  /(n  -  l)b!  +  b2]  -  F 


n=  1 


N-l 


=  \  [(2N  -  l)bj  +  2b2]rT  -  j  (bj  +  2b2)rA  -  b1  £  /nbj  +  b2  -  F, 


n=l 

,  ,  <rT  "  rA)  (rT  +  rA)  1  t.1  , _ 

"  2  I  3  (rT  "  rl>  - - V1 - ~  -  7  bl  E  Abl  +  b2 


[t 


n=  2 


and  for  the  equal-annuli  midpoint  algorithm  the  force  error  is 

N 

E 

n=l 
N 

E 

n=  1 


1 

EEM  =  4  £  [rn(+)  +  rn(-)]2[rn(+)  ~  rn(-)]  "  FC 

n=l 
1  N 

=  trn(+)  ~  rn(-)Hrn(+)  +  rn(_}]  -  Fc 


f  b!  ^v^bT 


N- 1 


+  b£  +  »4>2  +  2  /nb  i  +  \>2  1  “  F 


n=l 


)■ 


1  3  ,  (r|  -  r?)(rT  +  rA)  x  N"1  _ 

=  --j  (rT  -  rA)  +  - -  +  J  bl  E  /nbl  +  b2 


n=  1 


Hence , 


IEEmI  =  2  ^EEql 


(A8) 


(A9) 


(A10) 
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APPENDIX  B 


NUMERICAL  EXAMPLE 


With  five  segments,  the  resultant  forces  are  here  compared  by  using  the  sug¬ 
gested,  alternate  algorithm  (equal-annuli  midpoints)  and  the  extant  algorithm  (pure 
equal-annuli  spacing).  The  applicable  summation  algorithms  are  different,  as  noted 
in  the  text. 

For  five  segments  the  segment  extremities  (the  same  for  both  algorithms)  are 
given  from  equations  (87)  as 


rl (-)  =  rA  =  6-45 
ri(+)  =  r 2  (-)  =  14.632942 
r2(+)  =  r  3  (-)  =  19.663252 
r  3  (+)  =  ri,(_)  =  23.646648 
r4(+)  =  r5(_)  =  27.049666 
r5(+)  =  rT  =  30.07 


From  equation  (94),  for  equal-annuli  midpoint  spacing  the  observation  radii  are 


r*  =  10.541471' 
r2  =  17.148097 
r3  =  21.148097 
r4  =  25.34815 
rs  =  28.559833 


(B2) 


which  enable  us  to  compute  the  force  coefficients  from  equations  (62): 

a!  =  8.6617041 
a2  =  4.2898194 
a 3  =  4.10003 
a4  =  3.0973926 
a5  *  3.4710539, 


(B3) 


Thus,  for  equal-annuli  midpoint  spacing  with  five  segments  the  force-summation  algo¬ 
rithm  is  given  by: 
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Fgn  =  8. 6617041f(10. 541471)  +  4. 2898194f (17. 148097)  +  4. 10003f (21.65495) 

+  3. 0973926f (25.348157)  +  3.4710539f (28. 559833)  (B4) 

For  the  pure  equal-annuli  algorithm,  as  in  the  extant  model,  the  observation 
radii  are  given  from  equation  (79)  as 

r  i  =  11.307641' 

r2  =  17.331567 

r 3  =  21.746349  '  (B5) 

r4  =  25.4052 
r 5  =  28.599732, 

so  that  in  accordance  with  equation  (49),  using  the  differences  [rn(+)  -  rn(_)],  the 
extant  model  with  five  segments  produces  the  force-summation  algorithm: 

F£q  =  8. 182942f(ll. 307641)  +  5. 03031f (17.331567)  +  3. 983396f (21. 746349) 

+  3. 40 31858 f (25.4052)  +  3. 0203342f (28 . 599732)  (B6) 

For  a  force  of  order  k  in  radius  the  correct  integration  process  is 


Fc  =  /  T  rk  dr  =  kTT  (rT+1  "  rA+1; 


and,  therefore,  the  relative  force  errors  in  these  two  algorithms  are  given  by: 


= 


These  relative  errors  are  plotted  in  figure  22,  and  they  demonstrate  why  the  alter¬ 
nate  algorithm  delivers  less  error  for  segment  curve  characteristics  up  to  at  least 
cubic  order  (k  =  3)  in  radius. 
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TABLE  1.-  COMPUTED  rpm  WITH  VEHICLE  VELOCITY 


Vehicle  velocity,  knots 

rpm,  rad/sec 

0-180 

22.1416 

190 

21.9424 

200 

21.3972 

210 

20.8520 

220 

20.3069 

230 

19.7617 

240 

19.2165 

250 

18.6714 

260 

18.1262 

270 

17.5811 

280 

17.0359 

290 

16.4907 

300 

15.9456 

TABLE  2.-  UNIFORM-SEGMENT  EVALUATION  RADII 


Evaluation  point,  n 

Evaluation  radii  for 
number  of  evaluation  points,  N 

3 

5 

7 

10 

1 

0.1667 

0.1 

0.0714 

0.05 

2 

.5000 

.3 

.2143 

.15 

3 

.8333 

.5 

.3571 

.25 

4 

.7 

.5000 

.35 

5 

.9 

.6429 

.45 

6 

.  7857 

.55 

7 

.9286 

.65 

8 

.75 

9 

.85 

10 

.95 

TABLE  3.-  UNIFORM-SEGMENT  ALGORITHM  RELATIVE  ERRORS 


Number  of  evaluation  points,  N 

Relative  error  in  blade  moment 
if  the  force  profile  is: 

Constant 

Linear 

Quadratic 

1 

0 

-0.2500 

-0.5000 

2 

0 

-.0625 

-.1250 

3 

0 

-.0278 

-.0556 

4 

0 

-.0156 

-.0312 

5 

0 

-.0100 

-.0200 

6 

0 

-.0069 

-.0138 

7 

0 

-.0051 

-.0102 

8 

0 

-.0039 

-.0078 

9 

0 

-.0031 

-.0062 

10 

0 

-.0025 

-.0050 
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TABLE  4.-  EQUAL-ANNULI  EVALUATION  RADII 


•  — — 

Evaluation  point,  n 

Evaluation  radii  for 
number  of  evaluation  points,  N 

3 

5 

7 

10 

1 

0.4082 

0.3162 

0.2673 

0.2236 

2 

.7071 

.5477 

.4629 

.3873 

3 

.9129 

.  7071 

.5976 

.5000 

4 

.8367 

.7071 

.5916 

5 

.9478 

.8018 

.6708 

6 

.8864 

.7416 

7 

.9636 

.8062 

8 

.8660 

9 

.9220 

10 

.9747 

TABLE  5.-  EQUAL-ANNULI  ALGORITHM  RELATIVE  ERRORS 


Number  of  evaluation  points,  N 

Relative  error  in  blade  moment 
if  the  force  profile  is: 

Constant 

Linear 

Quadratic 

1 

0.4142 

0.5000 

0.4142 

2 

.2144 

.1893 

.1145 

3 

.1446 

.1061 

.0537 

4 

.1091 

.0702 

.0313 

5 

.0876 

.0508 

.0206 

6 

.0732 

.0390 

.0146 

7 

.0628 

.0311 

.0109 

8 

.0550 

.0256 

.0085 

9 

.0490 

.0216 

.0068 

10 

.0441 

.0185 

.0056 

TABLE  6.-  EQUAL-ANNULI  MIDPOINT  EVALUATION  RADII 


Evaluation  point,  n 

Evaluation  radii  for 
number  of  evaluation  points,  N 

3 

5 

7 

10 

1 

0.2887 

0.2236 

0.1890 

0.1581 

2 

.6969 

.5398 

.4562 

.3817 

3 

.9082 

.7035 

.5946 

.4975 

4 

.8345 

.7053 

.5901 

5 

.9472 

.8005 

.6698 

6 

.8855 

.7409 

7 

.9629 

.8056 

8 

.8655 

9 

.9216 

10 

.9743 

STARBOARD 


Figure  3.-  Standard  blade  orientation. 
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Figure  5.-  Lift-force  profiles. 
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Figure  5.-  Continued. 
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Figure  5.-  Continued. 
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Figure  5.-  Continued. 
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Figure  5.-  Continued. 
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Figure  5.-  Continued. 
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Figure  5.-  Concluded. 
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Figure  6.-  Continued. 
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Figure  6.-  Continued. 


60 


FORCE  PROFILE,  Ib/ft  FORCE  PROFILE,  Ib/ft 


O  BLADE  1  O  BLADE  4 
a  BLADE  2  *  BLADE  5 

A  BLADE  3 


50 


40 

30 

20 

10 

0 


-10 


-20 

4  8  12  16  20  24  28  32 

BLADE  RADIUS,  ft 


(q) 


75.0 

-  *  j  ♦  • 

’  "  ■  ■  ]  1  -T-  T  -  -  - 

.  .  j  .  ;  :  .  |_a.  . 

170  knots 

/ 

62.5 

. t- 

a 

£ 

50.0 

. 

uj 

_l 

LL 

37.5 

;  1 

o 

oc 

Q. 

UJ 

25.0 

o 

oc 

.... 

,  /  * 

,  .  •  1  \  * 

o 

LL 

12.5 

.  i  J 

*  »  »  * 

•  *  i  -  *  i. '■ 

*  *■-+*  .*  2  ~  ; 

Or 

> 

♦  :  *1*  ~ 

i- 

.  1-  "  .  ; 

.  .i  .  .  .  ;  -.4  . .  . 

-12.5 

l*  .  .  1 

1  4  .  i.  .  1 

4 

8  12 

16  20  24  28  32 

BLADE  RADIUS,  ft 
(r) 


1  i 

gg  - -  - -  4  t— L- 

.180  knots  .  .  i  i  .  ! 


75 - - 

i 


(s) 


(t) 


Figure  6.-  Continued. 
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Figure  6.-  Continued. 
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Figure  6.-  Concluded. 
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Figure  7.-  Aerodynamic  coefficient  logic. 
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Figure  8.-  Simplified  stall  onset  radii,  azimuth 
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Figure  9.-  Total  blade  aerodynamic  forces  and  moments  excluding  tip  loss. 
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Figure  9.-  Concluded. 
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Figure  10.-  Total  blade  aerodynamic  forces  and  moments  including  tip  loss. 
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Figure  10.-  Concluded. 
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Figure  11.-  Basis,  total  rotor  outputs,  body  axis  system. 
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Figure  12.-  Total  rotor  outputs  of  extant  model,  constrained  condition. 
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Figure  13.-  Aerodynamic  errors  of  extant  model,  constrained  condition, 

five  segments. 
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Figure  13.-  Concluded. 
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Figure  14.-  Total  rotor  outputs  of  extant  model,  balanced  condition. 
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Figure  16.-  Aerodynamic  errors  of  alternate  model,  constrained  condition 

five  segments. 
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Figure  16.-  Concluded. 
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Figure  17.-  Total  rotor  outputs  of  alternate  model,  constrained  condition. 
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Figure  18.-  Total  rotor  outputs  of  alternate  model,  balanced  condition. 
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Figure  19.-  Torque  comparison  of  alternate  model  with  five  evaluations. 
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Figure  20.-  Lagging  coefficient  variation 
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